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ABSTRACT 


The  purpose  of  this  thesis  is  to  obtain  the  approximate 
distribution  of  an  estimate  of  the  serial  correlation  coefficient 
of  a  prestationary  linear  Markov  process.  The  relevant  literature 
is  reviewed  in  Chapter  I.  Chapter  II  describes  the  prestationary 
linear  Markov  process,  gives  Daniels'  development  of  the  distributions 
of  a  sample  serial  correlation  coefficient  for  the  stationary  linear 
Markov  process  with  known  and  unknown  means  and  reviews  Patton's 
stationary  distribution  of  an  estimate  of  the  serial  correlation 
coefficient.  In  Chapter  III,  following  Daniels'  method,  we  obtain 
the  approximate  distribution  of  an  estimate  of  the  serial  correlation 
coefficient  when  the  mean  is  known  and  in  Chapter  IV  for  unknown  mean. 
We  also  compare  our  distributions  with  those  of  Daniels'  for  p  =  0. 
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CHAPTER  I 

INTRODUCTION  AND  SUMMARY  OF  LITERATURE 


In  the  present  work  our  object  is  to  find  an  approximate  distri¬ 
bution  of  a  serial  correlation  coefficient  in  a  pre-stationary  linear  Markov 
process.  Though  serial  correlation  has  attracted  many  persons  yet  no  one  has 
so  far  taken  up  the  problem  in  a  pre-stationary  case. 


Serial  correlation  tests  and  distributions  have  been  studied  in 
various  cases.  Below  we  give  a  summary  of  some  literature  on  the  subject. 

Anderson  [l]  defined  the  lag  L  serial  correlation  for  N  observa¬ 
tions  X^,  X^,  . . . ,  X^  by  the  expression 


lrn  ■ 


X1XL+1'1''  *  •'I'XN-lWl'(ZX1^^ 

ZXj  -  (EX.)2/fo 


Under  the  assumption  that  X^,  X^,  ...,  X^  are  independent  N(0,l)  he 
derived  the  exact  and  large  sample  distributions  of  ^R^  and  gave  a  table 
of  the  exact  significance  points  for  all  values  of  N  <  70.  For  N  >  70 
the  large  sample  approximation  could  be  used.  He  proved  that  the  distri¬ 
butions  of  and  ^R^  were  the  same  if  L  and  N  were  prime  to  each 

other  and  stated  that  the  distribution  of  ^R^  could  be  obtained  for  any 
L  and  N  if  it  were  known  for  values  of  L  and  N  for  which  L  was  a 

factor  of  N.  He  also  derived  the  distribution  of  ^R^  and  gave  the  sign- 

N 

ificance  points  in  the  cases  when  —  =  2,  3  anc^ 

Koopmans  [19]  considered  the  stochastic  process 


Xt  "  D  Xt-1  +  Zt 


(t=l,2,3, • • • ) 


•  ■  ie  -i  ..  .  .  .  .  ' 
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where  z were  independent  drawings  from  a  normal  distribution  with  mean  zero 
and  variance  cr2  and  the  coefficient  p  an  unknown  constant  such  that 
|p|  <  1.  On  the  basis  of  the  observed  values  x^,  x^,  ....  xT>  the  maximum 
likelihood  estimate  p  taken  as  an  estimate  of  p,  was  shown  to  be  a  function 
of  the  three  quadratic  forms 

o  2  2 

a  =  x^^  +  xT 

m  =  x^  +  . . .  + 

2  2 

and  n  =  x2  +  . . .  +  xT_^ 


For  testing  the  hypothesis  that  p  =  0  the  author  showed  that  it  was 

sufficient  to  know  the  distribution  of  and  in  this  case  studied  the 

l+n 

distribution  of  r  =  ^  where 

P 


2  2  2 
p  =  +  x2  +  . .  .  +  xT 


and  q  is  a  quadratic  form  in  the  variables  x^,  x2>  x^,  with  charact¬ 
eristic  values  k^,  k2>  k^.  Assuming  that  the  variates  x^,  x2>  x^ 

are  independent  N(0,l),  he  derived  an  expression  for  the  probability  density 
h(r)  of  r  taking  q  =  m.  He  also  considered  the  case 


q  =  ra  =  m  +  x^ 


and  showed  that  the  characteristic  values  of  m  and  m  were 


kt  -  frr 


r  2?rt 

kt  =  cos  — 


respectively.  He  obtained  an  approximate  formula  for  the  probability  density 

of  r  =  —  as 
P 


*3wa  3rre3anoo  mon'au 


-  3  - 


T 

4(T-1)22  jt"1 


r 

j 


arc  cos  r 


o 


(cosa- 


2 

s  in-|Tas  iriada  . 


Dixon  [6]  considered  distributions  related  to  the  quantities 


n 


1 


where  xn+^  =  >  the  x's  being  independent  N(a,cr2)  and  !>  being  the 

lag.  Using  the  relationship 


the 


x  = 
a 


likelihood  criterion 


a  +  bV< 

to  test 


(a) 


H, 


b  =  0  ,  was  found  to  be 


^  =  (l-b2)  ,  where 


A 

b 


A/v. 


n'  n 


and 


-  k  - 


(b) 


0 


,H1(a=0) :  b  =  0  ,  was  found  to  be 


n 


..  /  A  Ov2  A 

A  =  (l-b2)  ,  where  b  = 

01  '  0;  0 


E  x  x 
g  g+l 

E  x2 
g 


Approximations  to  the  moment  generating  function  were  used  to  derive  the 

/*  A 

moments  of  b  and  b^  which  were  found  to  be  exact  for  moments  less  than 

—  ,  where  g  is  the  greatest  common  factor  of  i  and  n.  These  moments 

g 


A  A  ~  2/n 


2/n 


were  used  to  derive  the  distribution  functions  of  b,  b^,  A^'  and  ^A^ 

A 

For  the  distribution  of  b,  a  Pearson  Type  I  approximation  was  used  to 
determine  Vj0  and  5$  positive  and  negative  significance  levels.  The  moments 
and  distribution  were  also  obtained  for 


and 


O'1! 


r  2 

a 


z  xf 

l 


2  = 


n 


and  the  first  2  moments  for 


2(l-b0) 


n 


n.  =  E  (x. 
'2  -y  v  i 


-  2x .  +  x.  )2/  v 

l+l  i+2 '  '  n 


Dixon  also  set  up  the  general  A-criterion  to  test 


H  :  b  . ,  b  b  =0  for  x  =  a  +  E  b.x  , 

r,m  m+1  m+2  r  g  1  1  a"^L 


and  gave  the  mean  and  variance  for  r  =  2  and  m=0,  1  (a=0,  and 
a  /  o).  He  also  indicated  the  set-up,  without  solution,  for  the  serial 


correlation  in  several  variables. 


- 
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Rubin  [28]  proved  the  equality  of  approximations  given  by  Koopmans 
[19]  and  Dixon  [6]  to  the  distribution  of  the  serial  correlation  coefficient 
r  ,  assuming  that  the  true  value  p  of  r  was  zero,  where 


Xt+i  =  *1  and  xt  are  indePendent  N(0,ct2)  . 

Madow  [22]  extended  the  results  of  Anderson  [l]  and  derived  the 
distribution  of  the  serial  correlation  coefficient,  using  the  circular 
definition,  for  the  population  value  p  /  0  . 

Leipnik  [20]  considered  the  circular  sequence  of  random  variables 
Xq,  x^,  x^  =  Xq  such  that  the  differences 


(  | p |  <  1 ,  t  =  1,2,...,T)  were  independent  and  normal  with  zero  means  and 
equal  variances  cr2  .  As  estimates  for  a2  and  p  the  expressions  ^ 
and  r  =  ^  respectively  were  used,  where 

p  =  x2  +  x|  +  . . .  + 

and  q  =  XjX2  +  . . .  +  x^  , 

By  smoothing  the  joint  characteristic  function  of  p  and  q,  he  deduced 
appropriate  expressions  for  the  frequency  functions  of  p  and  f  in  the 
case  of  an  arbitrary  p  .  For  r  he  obtained  the  frequency  function 


R  (r)  =  r^T+l)  (l-r2)2^T  l)  (i+p2-2pr)  2T 

0  r(^)  r(iT+i) 


■ 
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He  exhibited  the  graphs  of  this  distribution  for  T  =  20  and  various  values 
of  D.  He  also  obtained  general  expressions  for  the  approximate  moments  of 
p  and  r  and  explicitly  evaluated  the  mean  and  variance  of  these  variables. 

Anderson  and  Anderson  [2]  tested  the  independence  of  z t  ,  which 
is  the  normal  disturbance  in  a  seasonal  variation  of  a  time  series  xt  , 
against  the  hypothesis  of  linear  regression  of  the  type 

zt  =  p  Zt-L  +  Ut  (p  /  °) 

by  the  likelihood  ratio  criterion  which  is  shown  to  be  equivalent  to  the 
circular  coefficient  R  .  They  also  obtained  exact  and  approximate  distri- 

Li 

butions  of  R  and  tabulated  its  significance  points  for  L  =  1. 

Li 

Ogawara  [2k]  used  normal  regression  theory  to  estimate  the  para¬ 
meters  in  a  stationary  Gaussian  Simple  Markov  Process  under  the  assumption 
that  alternate  values  of  the  process  were  fixed,  set  up  the  confidence  limits 
for  the  autocorrelation  coefficient  and  indicated  the  theory  for  a  process  of 
order  h. 

Moran  [2j]  suggested  testing  serial  correlation  among  the  resid¬ 
uals  from  the  regression  line  to  make  sure  that  the  apparent  serial 

correlation  of  the  dependent  variable  may  not  be  entirely  due  to  serial 
correlation  of  the  independent  variable.  He  considered  the  first-order 
cyclic  serial  correlation  coefficient 

_  5  Vi+i 

“  T-»  9 

L  n  p 

r? 

where  e  .  =  e,  ,  and  derived  EfR,),  the  expected  value  of  R. ,  and 
n+ 1  1  v  1 '  1 

V(R^)  ,  the  variance  of  R^,  under  the  assumption  that  the  true  residuals 


1 


> 


are  independently  distributed  in  the  same  normal  distribution  and  stated  that 


Rx  -  E^) 

is  asymptotically  normally  distributed. 

Watson  and  Durbin  [32]  presented  an  exact  non-circular  test  of 
the  existence  of  serial  correlation  in  a  series  of  n  observations. 


Quenouille  [26]  derived  the  joint  distribution  of  circularly  defined 
serial  correlation  coefficients  on  the  assumptions  that  the  observations  were 
normally  and  independently  distributed.  The  results  could  be  extended  to 
the  case  where  the  observations  satisfied  the  relation 


aoxi  + 


alXi-l  + 


+  a  x .  =  e  . 

m  i-m  1 


where  e.  =  e.  and  the  e's  are  normally  and  independently  distributed, 
l+n  1 

Since  the  distributions  were  difficult  to  compute,  only  approximate  forms 
were  investigated. 


Gibson  [12],  in  an  investigation  which  involved  matching  a  series 
of  tree-ring  width  measurements,  constructed  a  routine  method  for  the  rapid 
measurement  of  the  lagged  correlations  between  two  series.  Instead  of  the 
usual  second  order  moment  correlation  he  used  a  fitst  order  moment  correl¬ 
ation  based  on  absolute  values  of  sums  and  differences. 


Durbin  and  Watson  [8]  taking 
Y  =  X  3  +  e 

where  Y  was  the  n-component  column  vector  of  dependent  variates,  X  the 
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(n  x  k)  vector  of  independent  (fixed  variates),  3  the  vector  of  regression 
coefficients  and  €  the  vector  of  (unobservable)  random  variables,  tested 
the  hypothesis  that  the  components  of  €  were  uncorrelated  against  the 
hypothesis  that  they  were  serially  correlated.  They  used  the  statistic 

z  'A  z 

r  =  “ 

where  A  was  some  given  positive  definite  matrix  and  z  =  M  y  was  the  vector 
of  residuals  of  an  observed  vector  y  from  the  fitted  regression  on  X  and 
M  =  I  -  X(X'X)  ^  X  .  The  ratio  of  quadratic  forms 

r=yf  M'  AMy  /  y1  MMy 

could  be  reduced  to 


by  an  orthogonal  transformation  y  =  H  £  .  If  the  components  of  e  were 
independently  N(0,cr2)  it  was  shown  that  £  ,  £  ,  . ..,  £  ,  were  also 

independently  N(0,cr2)  and  the  inequalities  on  the  distribution  of  r 
were  obtained  by  obtaining  the  inequalities  on  v^,  • • • »  vn  ^  *n  terms 

of  the  latent  roots  of  A.  In  [9]  upper  and  lower  bounds  to  the  5$>  2.5 $ 
and  1$  significance  points  of  the  distribution  function  of  the  criterion  in 
[8]  were  presented  for  1  to  5  fixed  variates  in  the  regression  model 
and  for  15  to  100  observations  and  an  approximate  procedure  of  testing  was 
also  outlined  when  the  bounds  test  was  inclusive. 

Sastry  [29]  gave  the  expressions  for  the  expected  value  of  the 
bias  in  the  empirical  autocorrelation  coefficient  caused  by  the  correction 
for  the  mean  and  obtained  some  numerical  results  in  some  representative 


cases . 


. 
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Walker 


•  •  •  I 


X 

n+u 


[31]  considered  n+u  (u  >  0)  consecutive  observations 
,  from  the  stationary  autoregressive  process 


i=o 


Y 

t 


where  {Y^}  is  a  stationary  m-dependent  process  with  E(Yt)  =  0 
E(Y^)  <  00  and  showed  that  the  distribution  function  of 

n 


Z 

n ,  u 


"'*1  (Xt  xt+u  -  E<xt  Xt+U» 

t=l 


and 


tends  to  that  of  a  normal  variate  with  zero  mean  and  finite  variance. 


Jenkins  [I7]  proposed  to  stabilise  the  variance  of  the  first  auto¬ 
correlation  coefficient  r  (of  a  sample  series  generated  by  a  first-order 
autoregression  with  E(r)  =  p)  by  the  use  of  sin  ^(r)  and  showed  that 
this  variate  was  asymptotically  normal,  more  stable  than  tanh  (r)  and 
useful  for  test  purposes  when  |p|  <r  0.9* 

Rao  [27]  considered  the  equation 

x  =  E(x  .)  +  €.,  1  <r  s  <  k,  1  <  t  <r  n 

st  v  st7  st  —  —  —  - 

where  E(x  )  is  a  linear  function  of  lagged  values  of  k  variables  in  a 
s  t 

simultaneous  equation  model  with  k  endogenous  variables  and  n  time 
periods  and  tested  the  hypothesis  that  the  n  normal  vectors  are 

independent  against  a  wide  class  of  alternatives,  where 

€t  =  f€lt’  €2t’  €st’  €kt^  * 
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Hannan  [14]  showed  that  Ogawara's  [24]  exact  test  was  asymptotically 
fully  efficient  if  the  parent  scheme  is  a  first  order  autoregression  and 
extended  this  test  to  give  an  exact  test  of  independence  of  residuals  from 
a  regression  equation  with  given  values  of  the  regression  coefficients. 

Watson  [33]  considered  the  consequences  of  assuming  an  incorrect 
covariance  matrix  for  the  residuals  from  regression  of  a  set  of  observations 
and  deduced  inequalities  on  the  bias  in  the  estimates  of  variance  of  the 
regression  coefficients,  on  the  efficiency  of  the  estimates  of  the  regression 
coefficients  and  on  the  significance  points  of  the  various  t-  and  F-  tests. 

Daniels  [5]  applied  his  adaptation  of  the  method  of  steepest 
descent  to  improve  upon  earlier  results  especially  with  regard  to  the  order 
of  approximation  and  to  treat  more  general  cases.  He  derived  distributions 
for  the  first  order  serial  coefficient  in  a  stationary  Gaussian  Markov  process 
in  the  circular  and  non-circular  cases  and  for  the  joint  distribution  of 
the  first  m  ordinary  and  partial  coefficients  in  circular  autoregressive 
process  of  order  m.  He  considered  both  the  cases  of  known  and  unknown 
means.  He  obtained  the  interesting  result  that,  with  a  Markov  Process, 
Leipnik’s  [20]  formula  for  the  circular  case  needs  only  a  change  in  the  degrees 
of  freedom  to  apply  in  the  non-circular  case,  with  an  error  of  order  0(T 
where  T  is  the  number  of  observations. 


Watson  [34],  using  the  characteristic  function  and  latent  root 

tKL 

representation  of  the  serial  coefficients,  derived^approximate  joint  distri- 


CL. 


bution  for  the  case  of^/Gaussian  process  with  zero  autocorrelations  and  non¬ 
zero  means  and  compared  this  with  more  precise  results  by  Jenkins  [17]  and 
Daniels  [5]. 


■  :  •  '  -  •  ri  :  .  .9 
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Hannan  [15]  extended  the  tests  for  Serial  Correlation  [14]  to 
vector  Markov  Processes  and  tested  the  independence  of  successive  vector 


variates  by  examining  the  canonical  correlation  between  and 

He  also  extended  the  method  to  provide  an  exact  test  for 


the  serial  correlation  of  a  vector  or  residuals  from  regression  and  obtained 
bounds  for  the  asymptotic  efficiencies  of  the  tests.  He  [16]  also  obtained 
the  asymptotic  distribution  of  a  statistic  r  to  be  used  in  testing  for 
serial  correlation  in  the  residuals  of  a  regression  of  the  form  Y  =  X  (3  +  e 
with  the  components  of  e  generated  by  a  stationary  process.  He  also  found 
approximate  lower  bounds  for  percentage  points  of  the  distribution  of  r 
when  the  regression  was  on  orthogonal  polynomials  or  on  orthogonal  polynomials 
and  random  variables. 


White  [37]  considered  the  distribution 


_  N  N-l 


f(t)  =  C  (l-2pt+p2)  2  (1-t2)  2  (-1  <  t  £  1) 


which  is  an  approximation  to  that  of  the  serial-correlation  coefficient  of 
a  circular  first  order  Gaussian  autoregressive  process  and  derived  formulas 
for  its  moments. 


Durbin  [7]  extended  the  results  of  Durbin  and  Watson  [8  and  9]  t0 


models  with  p  +  1  interrelated  equations  in  p  +  1  dependent  or  endogen¬ 
ous  variables  (y)  and  k  +  q  (q  >  p  and  k  >  O)  independent  or 
exogenous  variables  (x).  Taking  simultaneous  equations  of  the  form 


Ay  =  B  x  +  € 


where  €  represented  the  vector  of  errors  or  disturbances  and  assuming  that 
n  simultaneous  observations  were  obtained  on  the  variables  and  that  q  of 


' 


,b 3930*0  y*mj  *a  c  yd  fears* on  9  So  :■  nai  ,«jtoo  a  3  dJj  v 


12 


the  elements  of  the  corresponding  vector  of  B  were  zero,  results  to  test 
for  serial  dependence  in  a  given  equation  of  the  system  were  obtained. 

r  1  11 

Leipnik  [21 J  used  a  method  of  Kac  and  Kac  and  Erdos  to  obtain  a 
closed  form  for  the  joint  generating  function  of 


where  the  X^'s  are  serially  correlated  normal  variables.  He  also  computed 

the  exact  moments  of  several  estimates  of  the  variance  and  autocorrelation 
Ifct 

in  the  case  of j^Ornstein-Uhlenbeck  process. 

Siddiqui  [30]  showed  that  for  the  first  order  autocorrelation 
coefficient  r  of  N  independent  N(0,l)  variables,  the  expansion  of 
Prob  (r  >  r^)  in  powers  of  1  -  r  begins  with 

(N-2) 

C(1  -  r  )  2  ,  O.lA  c  C  <r  0.22  . 

Weinstein  [36]  tested  various  estimators  for  the  serial  correla¬ 
tions  of  small  lags  in  a  stationary  series  with  an  unknown  mean  and  variance 
using  an  artificially  constructed  second  order  autoregressive  process  with 
samples  of  14  and  showed  that  the  result  of  eliminating  the  unknown  mean 
was  more  crucial  than  eliminating  the  unknown  variance. 

White  [38]  showed  that  the  limiting  distribution  (suitably 
normalised)  of  the  least  squares  estimate  a  of  the  autoregression  was  Cauchy 
when  |  a  |  >1  >  normal  when  |a|  <  1  and  of  a  third  form  whose  char¬ 

acteristic  function  was  obtained  by  him  when  |a|  =  1.  He  [39]  further 
showed  that  for  any  a,  with  the  possible  exception  of  the  case  | a |  =  1, 


the  statistic 


■ 
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is  asymptotically  N(0,l)  where  a  and  u  are  the  maximum  likelihood 
estimates  of  the  parameters  of  the  autoregression 


=  a 


Xt-1  + 


u 


the  u's  being  independent  N(0,ct2).  He  [40j  also  obtained  expansions  for 

A 

the  expectation  and  variance  of  the  maximum  likelihood  estimate,  a,  of 

the  serial  correlation  coefficient,  a,  of  a  first  order  autoregressive 

-3  b  , 

Gaussian  process,  as  far  as  the  terms  in  T  and  a  (where  T  is  the 
length  of  the  series)  for  the  case  where  the  first  observation  was  uncondi¬ 
tional  and  was  known  to  be  equal  to  zero. 


Griliches  [13]  considered  the  model 


y„  =  3  x+ 


7  7t-i  + 


=  °et-i + 


w 


where  w  are  independently  distributed  random  variables  with  zero  means 
and  y^  and  e^  are  stationary  stochastic  processes  and  showed  that  the 
simple  least  squares  estimator  c  of  7  was  biased  and  that  the  bias 
was  positive.  He  further  showed  that  if  7=0  but  the  estimating  equation 


t  =  b  xt  +  c  yt_1 


was  used  to  reduce  the  serial  correlation  of  the  residuals,  the  expected 


0(1 


) 


value  of  c  would  be 


14  - 


Chanda  [3]  proved  that  the  maximum  of  the  absolute  value  of  the 
sth  order  serial  correlation 


was 


and 


n-s 


n 

s  -  n-s 


x*. 

t  t+s 


n 


y 


n  tc 

-  cos  — - 

n-s  m+2 


for  n  =  ms  +  u  ,  u>0 


m  it 

— 7  cos  — 7 
m- 1  m+ 1 


for  n  =  ms 


In  the  present  thesis  we  consider  the  linear  Markov  process 


s-1  s-2 

x  =p  e.  +  p  e_  +  . . .  +  p  e  ,  +  e 

s  1  2  s- i  s 


where,  e  ,  e  ,  e  are  independent  N(0,l)  random  variables  and  p 

I  d  s 

(the  serial  correlation  coefficient)  is  such  that  |p|  <"  1.  We  review 
Daniels'  method  and  work  in  obtaining  the  approximate  distribution  of  the 
sample  serial  correlation  coefficient 


where 


c  =  X,X_  +  .  .  ,  +  X  ,x 
1  2  n-1  n 


C0  =  M  +  X2  +  +  Xn-1 


and  then  use  Daniels'  method  to  obtain  the  approximate  distribution  of  r 
in  the  prestationary  linear  Markov  process  under  consideration.  We  obtain 
the  distribution  of  r  in  the  case  of  known  as  well  as  the  unknown  mean. 


:  “J  '  • f.  •: 


' 


CHAPTER  II 


THE  APPROXIMATE  DISTRIBUTIONS  OF  ESTIMATES  OF  THE  SERIAL 
CORRELATION  COEFFICIENT  FOR  A  NON -CIRCULAR  MARKOV 
PROCESS  WITH  KNOWN  AND  UNKNOWN  MEANS 


1 .  Linear  Markov  Processes 


Consider  a  subset  of  time  points 


t 


1 


c  t2< 


...  <  t 

n 


and  the  random  variables  X(t.)  at  time  t.. 

J  J 

X( t^) ,  ...,  X(t^)  is  denoted  by  p[X(t1),  ... 
stationary  process  if 


The  joint  distribution  of 
,  X(tn)].  We  have  a  completely 


p[X(t1),...,X(tn)]  =  p[X(t1-T),...,X(tn-T)] 

for  all  subsets  t^  <  t^  <  ...  <  t^.  This  means  that  the  joint  distribution 
is  invariant  under  a  change  of  time  origin.  A  process  is  stationary  of 
order  k  if 

£.  [xa'1(t1)  ...  xan(tn)]  =  gtx’V^r)  ...  X  n(tn-r)] 

for  all  t  and  all  a,  +  a_  +  . . .  +  a  <  k  . 

12  n  — 

We  shall  now  consider  the  cases  of  order  one  and  two.  If  the 
process  is  stationary  of  order  one, 


£[X(t)]  =  0  [X(t-r)  ]  =  u 

for  all  r  which  shows  that  the  mean  is  independent  of  time.  If  the 


process  is  stationary  of  order  two,  we  have 


for  all  t  ,  and 
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£[X2(t)]  =  £[X2(t-r)) 

gWtjJXftg)]  =  gW^-rWlyr)] 

for  all  r  .  Then 

Var  [X(  t)  ]  =  £[X2(t)  ]  -  {£  [X(  t )  ] }  2 

=  £[x2(t-T)M£[x(t-T)]}2 

=  Var [X( t--r) ] 

=  <T2 

for  all  t  .  This  shows  that  the  variance  is  independent  of  time.  We  call 
cr2  the  autovariance.  Also 

cov[x(tl),x(t2)]  =  £[x(tL)x(t2)]  - £[x(t1)]  g[x(t2)] 

*  £[X(t1-T)X(t2-T)]  -  £[X(t1-r)]  £[x(  t2«r)  ] 

=  Cov[X(t1-r),  X( tg-r ) ] 

- 

for  all  t  .  C(t2-t^)  is  called  the  autocovariance  and,  as  shown,  is 

independent  of  the  time  origin. 


The  autocorrelation  function  p(t')  is  defined  as; 

C<V 

cr2 

where  r'  =  t  -  t^  is  the  lag.  p(r')  measures  the  degree  to  which 
variables  at  time  difference  r'  are  correlated,  p(t')  is  an  even 
function  of  t' .  It  can  also  be  shown  that  necessary  and  sufficient 


conditions  for  a  function  to  be  an  autocorrelation  function  are  that 
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p(o)=l,  p(r')  is  even 

ant^  p(r')  is  positive  semi-definite, 

A  process  which  is  stationary  of  order  one  but  not  stationary  of 
order  two  will  be  called  a  pre-stationary  process. 


The  linear  Markov  process  in  discrete  time  is  defined  by 
(2.1.1)  xs  =  D  xs_i  +  es  (s  =  1,2, ...,n) 

where  the  constant  p  is  called  the  serial  correlation  coefficient  and  the 

{e  )  are  independently  and  identically  distributed,  N(0,l)  random  variables, 
s 

Following  Patton  [  (25)  ,pageS9"H]  we  have 


s 

P 


X  — ►  o 

o 


where  x  is  the  initial  state, 
o 


and 


Var (x  )  = 
v  s' 


y 


COV(XSlXS-T) 


for  all  t  if  | p |  <  1  and  s  -*•<»  . 

Hence  if  |p|  <  1  ,  the  process  ultimately  becomes  stationary  to 
the  second  order,  which  in  the  case  of  a  normal  process  implies  complete 
stationarity . 

Now  let  us  consider  a  linear  Markov  process  in  discrete  time 
defining  xi  (i  =  l,2,...,n)  as  follows: 
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(2.1.2) 


X1  =  el 


X2  =  DG1  +  e2 


4  x  =  p2e  +  pe  +  e 

3  123 


n-1  n-2 

x=p  e.,  +  p  e^+.,,+pe  ,+e 

n  1  2  n-1  n 


where  e^  (i  =  l,2,...,n)  are  independent  N(0,l)  random  variables  and 

|p|  <  1  . 

We  observe  that  for  s  =  1,2,..., n 

s- 1  s-2  2 

x  =p  e.  +  p  e_  +  , . .  +  p  e  _  +  pe  .  +  e 

s  1  2  s-2  s-1  s 

Since 

£(ei)  =  0  ,  (i  =  1,2,  ...,n)  > 

£ (xg)  =0  ( s  =  1,2, ... ,n)  , 

Hence  the  process  is  stationary  to  the  first  order. 

Again  since 

£  (e2)  =  1  (i  =  1,2, . . . ,n)  , 

and 

£ (e±eJ)  =0  (i  /  j) , 


we  have 


2(s-l)  2( s-2) 

P  '  +  P  + 


4  2 

•  •  •  +  p  +  p  + 


1 


-  19  - 


and 


6<*J.T>  - 


1  -  P 


2s -2t 


1  -  D‘ 


Therefore 


£<x?)  /  £<xs-t)  • 


Again 


and 


£,  x  s+t-2  s+t-L  s-t 

(xs  xt)  =  P  +  P  +  ...  +  P 


(s  >  t) 


s-t  f  2t-2  _2t-l+  2  .  \ 

=  P  (P  +  P  +  . . .  +  p2  +  1) 


.s-t  f  1  -  p 

3  0  <  1^75 


2t 


£  <VT  Xt-T)  3 


s-t 


1  -  P 


2t-2T 


1  -  p‘ 


Hence 


£<XS  xt>  /  £(xs-t  Xc-T)  • 


Therefore  the  process  defined  in  equations  (2.1.2)  is  a  pre-stationary  one. 

2.  Non-circular  Markov  Process;  Daniels1  Problem, 

Consider  the  process  x  =  p  x  .  +  e  for  all  s  where  the  e's 

s  s-1  s 

are  independent  N(0,l)  random  variables.  We  assume  that  |p|  <  1  and  that 
the  process  is  stationary. 

Let  x, ,  x^,  ....  x  be  a  sample  of  observations  from  the  station- 
12  n 

ary  process.  They  have  a  joint  multivariate  normal  distribution  since 
e^,  e^,  ...,  en  are  independent,  N(0,l)  random  variables.  The  variance 


n 


covariance  matrix  of  x^,  x^,  ..., 


is  the  (nxn)  matrix 


20 


(2.2.1)  C  = 


- 

1 

P 

n-1 

p 

1  -  P2 

1  -  P2 

•  •  « 

1  -  P2 

p 

1 

n-2 

P 

1  -  p2 

1  -  P2 

•  •  • 

1  -  P2 

• 

D 

3 

i 

»  • 

• 

• 

o"-2 

1 

i  -  p2 

1  -  p2 

1  -  P2 

L 

n 


Multiplying  the  second  column  by  p  and  subtracting  it  from  the 
first  column  in  |c|  and  repeating  the  same  procedure  for  third  and  second 
columns  and  so  forth,  we  get 


(2.2.2) 


£i  ■  rh*  ■ 


-1 


It  is  easily  shown  that  C  is  the  (nxn)  matrix 


(2.2.3)  £ 


-1 


-P 


1+P‘ 


-P 


(0) 


-p 


l+p‘ 


(0) 


-p 


l+pe 


-p 


The  joint  distribution  of  x^,  x^,  . . . ,  x^  is  given  by 


■P 

1 


(2.2.1+)  dF(x1,x2> . . .  ,xn)  = 


(2tf) 


“X  exP{"2  x  C  ^  x)  dx^.  •  *dxn  > 
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where  C  is  the  variance-covariance  matrix  and 


s'  =  <xr  x2 . xn)  • 


For  the  case  under  consideration,  where  |C|  and  C  ^  are  given  by  equations 
(2.2.2)  and  (2.2.3)  respectively,  we  have 


dF 


1. 

(x1,...,xn)  =  '  exp{-^[x^+(l+D;?)(x|+.  .  .+x^_1)+x^-2D(x1xp+. .  .+xn_1x^)  ]) 


(2«) 


2 


dx, . . . dx 
1  n 


The  sample  estimate  of  p  is  taken  to  be 


(2.2.5) 

where 

(2.2.6) 

and 

(2.2.7) 


r  = 


c  =  x.x„  +  x  x_  +  ...  +  x  ,x  , 
12  23  n-1  n 


1  2  2 

c  =  2x7  +  xr  +  . 

o^l  2 


2  12 

+  +  2'X  . 

n-1  n 


The  joint  Moment  Generating  Function  of  c  and  c^  is 

T  c  +  Tc 

M(Tq,T)  =  £(e  °  °  ) 


(l-P2)2 


* co  T  c  +Tc-^x'C  ^x 


o  o 


00 


(2it)‘ 


dx. dx  , . .dx 
12  n 
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Now 


Vo+Tc i'S's  -  -itd-T0)xf  +  (i+p2-8To)(x|+...+x^1) 


+(1-to)x^-2(p+t)(XiX2+. . .+Xn_lXn) ] 


=  -jr  x'  S  x 


where  B  is  the  (nxn)  matrix 


(2.2.8)  B  = 


1-T 


(d+t) 


-(p+t)  1+p2-2T  -(p+T) 


(0) 


•(p+T)  l+p2-2T  -(o+T) 


(0) 


-(p+T)  l+pa-2To  -(p+T) 


- (p ! t)  1-T 


Then  M(Tq,T)  may  be  written  as 


M(T  ,T)  =  f  ...  T  e”2  -  dx  . . .  dx  . 

o  n  f  /  1  n 


-oo  ' 1  -oo 


(at)' 


There  is  an  orthogonal  transformation  that  will  diagonalise  B.  Making  this 
transformation  and  performing  the  integration  we  get 


(2.2.9) 


M(T  ,T)  = 


B 


7 


where  |b|  2  is  the  Jacobian  of  the  transformation. 


From  equation  (2.2.8)  and  using  (ll,20)  and  (11.22)  we  observe 


that 
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B  =  f 


(p-p)f  -  q5  an-2  ,  (g-p)£-i-  -r  ftn-2 

a  -  3  a  +  a  -  3  1 


-  q 


Lrz&JL:.^.  an"3  Ca-p)f  +  q5  Rn-3 
a  -  3  a  -  3  P 


where 


(2.2.10) 


a  +  3  =  P  » 
a3  =  q2  , 


<  p  =  1  +  o2  -  2T  , 


q  =  -  (p  +  T)  and 

f  =  (1  -  T  )  . 

'  o' 


Substituting  for  p  and  q  in  terms  of  a  and  (3  we  get 


tel  =  f 


[qf  -  a3  n- 

a  "  3  a 


2  flP.  _i  M  B 
a  -  3 


-] 


[ 


.  af  -  aB  n-3  a3  -  Bf 

-  a3  - - *  a  ^  ^ 3 

!  a  -  3  a  -  3 


or 


n-1 


(2.2.11)  |B|  =  (f  -  3)2  “ 


a  -  3 


(a  -  ()s  ? 


n-1 


a  -  3 


Let 


(2.2. 12) 


z  + 


,  l  +  02  -  2T 
_1 _ o  _  2 

z  P  +  T  '  q 


which  gives  us 


(2.2.13) 


z2  +  -^z+l  =  0 

q 
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The  roots  of  this  equation  are 


z  = 


-  E  + 

q 


l 

q 


P  -  N TP  -  4^  1 


and 


1 

z 


l. 

2 


_  £ 

q 


£f 

q2 


l 

q 

a 


P  +  V~p2  -  Iq2 

2 


or 


a 


£  _  P  +  T 
z  z 


(2.2.14)  <  and 


(3  =  “  qz  =  (d  +  T); 


Substituting  for  p,  q,  a,  0  from  equations  (2.2.10)  and  (2.2.14)  in 
equation  (2.2.11)  we  get 

(2.2.15)  |B|  =  (  1"T0"DZ+Tz)2  (p+l)^1.^) 


,p+T  ,  _  vp  .  _Nn- 1  n- 1  _ z _ 

*(  z  '  1  +  To)  (0+T)  Z  (o+T)(l-zS)  • 


Since 


+  P2  -  2T 


P  + 


z  + 


1  z2  +  1 


(l  +  p2  -  2T  )z  = 
v  o' 


(p  +  T ) ( 1  +  z2) 
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which  gives 


(1  -  Tq)z  +  (os  -  To)z  =  p  +  T  +  z2(o  +  T) 


which  in  turn  gives  us 


1  -  T 


(o  ±  T) 

z 


z(p  +  T)  -  (p2  -  T  ) 


and 


1  -  T  -  z(p  +  T)  =  -  (p2  -  T  )  . 

o  z  o' 


Substituting  for  these  values  in  equation  (2.2.15)  we  get 

l2  _.m  n-2 


B 


P  +  T 


/  -  2  rr.  \  /P+T  V  “  _ 1 

o'  z  '  (1-z2) 


z(p+T)  -  (p2-Tq) 


(p+t) 


n-2  z 


n 


"fl-z"5) 


or 


(2.2.16) 


“(l-z2) 


rf  2(P*-T0)l2 

r  p2-Tof  2„- 

L 1  p  + T I 

r  '  (p-t)  j  -  _ 

As  mentioned  by  Daniels'  [53  if  1-r2  is  not  small,  leaving  out 
2n 

the  term  in  z  within  bracket  in  equation  (2.2.16)  incurs  an  exponent¬ 
ially  small  error  in  the  final  approximation.  Hence 

2 


(2.2.17) 


B 


,  (p+T)n 

z(p2-Tq) 

“  zn(l-22) 

1  ‘  (P+T) 

Now  let 

(2.2.18)  u  =  Tq  +  r  T 
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so  that 


1 


z  + 


1  +  P2  -  2T 
_ o 

P  +  T 


1  +  P2  -  2u  +  2rT 

P  +  T 


or 


(2.2.19) 


1 

z  +  — 


z 


1  +  P2  -  2Pi  -  2u 
P  +  T 


Hence 


and 


(2.2.20) 


which  gives 


(2.2.21) 


and 


(2.2.22) 


1  -  2rz  +  z2 
z 


1  +  P2  -  2pr  -  2u 

P  +  T 


z(l  +  p2  -  2pr  -  2u) 

1  -  2rz  +  z2 


/  ( 1  -  Pz)(p  -  z)  +  2uz  \ 
[  1  -  2rz  +  z2  J 


T  _  (1  +  z2) ( u  +  pr)  -  rz(l  +  p2) 

o  1  -  2rz  +  z2 


By  use  of  equations  (2.2.18),  (2.2.20)  and  (2.2.21)  we  obtain 


1 


z<p2  -  Tq) 
P  +  T 


1  .  (1  -  2rz  +  z2) ( p2  -  u  +  rT) 
( 1  -  2pr  +  p2  -  2u) 


( l-2pr+pg-2u)  -  (p2-u) ( l-2rz+z2)  +  r f ( 1-pz) (p-z)+2uz ] 

1  -  2pr  +  p2  -  2u 


or 


(2.2.23) 


z(p2  - 

1 - 

P  + 


(1  -  pz)n  -  pr  -  (r  -  p)z]  -  u(l  -  z2) 
1  -  2Pr  +  P2  -  2u 
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Finally  equations  (2.2. 20 )  and  (2.2.23)  together  with  equation  (2.2. I7) 
give  us 

I  _  |  zn(  l-2Pr+P2~2u)n  f  ( 1-pz)  f  l-pr-(r-p)z]-u(  1-z2)  \2 

"  "  Zn(l-^)(l-2rz+z2)n  >  1  -  2or  +  0*  -  2u  f 

=  .2pr±P— — -  ((l-pz)[l-pr-(r-p)z]  -  u(l-z2)}2 

( 1-z2)  ( l-2rz+z2) 


Substituting  for  |b|  in  (2,2.9)  we  get 


(2.2.24)  M(u-rT,T)  ^ 


( l-p2)^( l-z2)^( l-2r  z+z2) 


2-1 

>2 


( l-2pr+p2-2u)  (( 1-pz)  [  1-pr -(r  -p)z ]-u(  1-z2) } 


(2.2.25) 


Differentiating  equation  (2.2.20)  with  respect  to  z,  we  get 

St  _  ( 1-z2)  ( l-2pr+p2-2u) 

SI  "  ( 1-2  rz+z5)2 


Combining  equations  (2.2.24)  and  (2.2.25)  we  have 


M(u-rX,T)  ||  „ 


f  l-pa)a(  l-2rz+za) 


2 


n 


( l-2pr+P2-2u;  (( 1-pz) [ 1-p r-( r-p)z ]-u(  1-z2) } 


Hence 


a 

5; 


[M(u-rT,T) 


n  _2 

_ (  1-z2)^2(  l-p2)2(  l-2tz+z2)^ _ 

2-1 

( l-2pr+p2-2u)2  { ( 1-pz) [ l-pr-(r -p)z ]-u( 1-z2) ) 


X 


n  -  4  + 


(l-z2)(  l-2pr+pg-2u) 

((1-pz) [1-pr- (r-p)z]-u( 1-z2)] 


} 


' 
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so  that 


a 


[M(u-rT,T)  g2] 


u=o 


a  -2 

n(  l-z2)^/^  l-p2)2(  l-2rz+z2)2 

~  a-i 

(l-2pr+p2)2  ( 1-pz) [l-pr-( r-p)z] 


* 


( l-z2)(  1-2PE+P2) 

( 1  -P  z)(l-pr-(r  -p)  z  ] 


Thus  equation  (l.2)  gives  us 


(2.2.26) 


h(r) 


—  n(l-PS)2  

n 

2jti(  l-2Pr+P2)2 


1 


n 


J 


P  <p(  z)  ( l-2rz+z2)2 


-2 


where 


(2.2.27)  <P(z) 


(l-zs)l/2 

( 1-pz) [ l-pr-(r-p)z] 


(l-z2)(l-2pr+pg) 

( 1-pz) [ l-pr-(r-p)z] 


The  subsequent  discussion  in  this  section  is  essentially  that 
given  by  Patton  [25]. 

1  1  x  -  2rT 

With  u  =  0  ,  the  mapping  z  +  —  =  - - - 

z  p+i 

is  easily  seen  to  consist  of  the  following  elementary  mappings: 


(2.2.28) 


A: 

B: 

C: 

D: 

E: 


V  = 
U  = 

s  = 


P  +  T 


V 


1  -  2pr  +  o‘ 


1_ 

2U 


Q  =  S  +  r 


2  +  I=  2Q 


dz 
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A  maps  the  T-plane  cut  along  the  parts  of  the  real  axis  exterior  to  the 

interval  -j"  ^  onto  the  V-plane  cut  along  the  parts 

f  -( l-2pr+p2)  l-2pr+p2  \ 

{  2(l+r)  *  2(l-r)  [  ' 


of  the  real  axis  exterior  to  the  interval 
B  maps  the  region  of  the  V-plane,  as  above,  onto  the  U-plane  cut  along 

maps 


f  - 1  1  1 

the  real  axis  exterior  to  the  interval  1  2( 1+r)  *  2( 1-r)  j  *  C 


this  region  of  the  U-plane  onto  the  S-plane  cut  along  the  real  axis  from 
-(1+r)  to  (1-r).  This  region  of  the  S-plane  is  mapped  onto  the  Q-plane 
cut  along  the  real  axis  from  -1  to  1  by  D,  and  this  region  of  the 
Q  plane  is  mapped  onto  the  interior  of  the  unit  circle,  |z|  =  1  ,  in 
the  z  plane  by  E.  Hence  the  net  effect  of  the  transformation  is  to 
map  the  whole  T-plane  cut  along  the  real  axis  exterior  to  the  interval 


fdLkelf  U-p)2 

l  2(l+r)  ’  2( 


-J>V  I 

r+y  f 


onto  the  interior  of  the  unit  circle. 


Since  by  equation  (2.2,19) 


1 

z  +  - 

z 


u=o 


1  +  Pg  -  2pr 

P  +  T 


+  2r  , 


as  T  -v  +  i  oo  ,  z  +  -  -■►  2r  . 
~  z 


If 


then 


and 


z  +  -  =  2r 
z 


z2  -  2rz  +1=0 


z  =  r 


+  i  sfrr 


+ie 

=  e— 


' 
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where 


r  =  cos  9 


We  also  note  that  if  T 


z  + 


=  0 

1  _ 

z 


and  hence,  the  transformed  path  in  the  z-plane  cuts  the  real  axis  at  z  =  p. 
Thus,  the  path  of  integration  in  the  z-plane  crosses  the  real  axis  at  z  =  p 

X0  iQ 

and  ends  on  the  boundary  of  the  unit  circle  at  e  and  e  where 
r  =  cos  9  . 


The  only  possible  singularity  in  the  integrand  of  the  integral  for 
h(r)  would  be  the  real  singularity  arising  from 

1  -  pr  -  (r-p)z  =  0  . 

This  singularity  is  avoided  if,  on  the  real  axis,  we  confine  z 
to  the  following  range; 

(i)  P  <  z  <  r  ,  P<r 

Now,  we  may  write 

1  -  pr  -  (r-p)z  =  (l-z)(r-p)  +  (l-r)(l+p) 

Then,  since  |p|  <  1  and  |r|  <  1  , 

we  see  that 

1  -  z  >  0,  r  -  p  >  0,  1  -  r  >  0  and  1  +  p  >  0  . 

From  this  it  can  be  seen  that 


(l-z)(l-p)  >  0 


and 


( l-r) ( 1+p)  >  0  , 


- 
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and  finally  that 


1  -  pr  -  (r-p)z  >  0 


(ii)  r  <_  z  <  o  ,  r  <  p 


Here, 


1  -  pr  -  (r-p)z  =  1-r2  +  (p-r)(z-r)  . 


As 


| P |  <  1  ,  and  | r |  <  1  , 


we  have 


1  -  r2  >  0  ,  p-r>0  and  z  -  r  >  0 


Then 

(p-r)(z-r)  >  0 

and  hence 

1  -  pr  -  (r-p)z  >  0 


( iii)  z  =  r  =  p 

Since  |p|  <  1  and  |r|  <  1  , 

1  -  pr  >  0  and  (r-p)z  =  0 


so  that 


1  -  pr  -  (r-p)z  >  0 


Since  there  is  no  singularity  in  the  closed  interval  joining  p  and  r, 
we  may  deform  the  path  of  integration  in  the  z-plane  to  be  the  straight 
line  joining  e  ^  to  e^  and  crossing  the  real  axis  at  z  =  r  . 


-  32  - 


On  the  path  of  integration  we  may  write 

X. 

(2.2.29)  z  =  r  +  i  a)  ( 1-r2)2 
where  u>  is  the  real  variable  with 

-1  <  u>  <  1  . 

Then 

(2.2.30)  1  -  2rz  +  z2  =  ( l-r2)(  1-u)2) 

and 


(2.2.31) 


dz  -  i(l-ra)2  dw  . 


Substituting  equation  (2.2.29)  in  equation  (2.2.26),  we  have 

n-3 

(2.2.32)  h(r)  (l'r2) -  f  (1-m2)^  du 

_ J-  1  J-1 


2«(  l-2pr+p2)' 


The  integrand  of  h(r)  is  of  the  form 


<p(z)  ty(z)]2 


where 


^(z)  =  1  -  2rz  +  z2  . 

Since  the  integral  cannot  be  readily  evaluated  in  closed  form,  we  expand 
Cp(z)  as  a  power  series  in  z  -  r  ,  where 


ip'(r)  =  0 


,C  +  - 


j  ■  .  '  -1  &r;  ' 


l  '  -v  +  • 


2. S)  $ak3una%4ot 


So  tti n^9ial  oifT 
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and  integrate  with  respect  to  w  .  Using  equations  (2,2.29)  we  have 

z  -  r  =  i  u  (l-r2)i 


and  hence, 


(k),  \  /  \k 

q>(z)  =  «p(r)  +  cp'(r)(z-r)  +  ...  +  2 -  k  +  ••• 


or 


.k  k. 


(2.2.33)  9(z)  = 


k 
-2\ 2 


cp(r)  +  cp'(r)  i  w(l-r2)^  +  ...  +  1--‘ ^  q/k)(r)  +  ... 


Substituting  this  equation  in  (2.2.32),  we  get 


n- 


(2.2.3k)  h(r)  ^  SIIl£f)lUzlf) - 

2n( l-2pr+p2)2 


00 

I 

k=o 


ik(l-r2)2  ®(k)(r)  T1  k,.  ^ 

— 1 - ^ ^  /  0)  (1-w2) 

'-1 


-  -2 

2 


da) 


Now,  if 


k  =  2m  +  1  ,  m  =  0,1,2,.,. 


then 


2 


wk  (1-w2)2 


is  an  odd  function  and  the  integral  vanishes.  Therefore,  equation  (2. 2.3)1) 
reduces  to 


n-3  w 

n( 1-Pa)^( l-rg)  2 


2rt( l-2pr+p2) 


dL±  y 

b.  1  L 

2\  2  k=o 


(-l)k(l-r2)k  <p(gk)(r) 

(2k)! 


B_  2 


J 


r  2k, .  2v2  , 

I  w  (l-w^)  dw  . 


-1 


(2.2.35)  h(r) 


t90tl  1  \ 


2 


-  3b  - 


To  evaluate 


(2.2.36) 


\  - 


f1  2k  .  2. 

=  /  w  (l-W^) 


i-2 


du> 


-1 


we  proceed  in  the  following  manner!  Taking 


u>2  =  u,  du)  =  \  u  ?  du  , 


we  have 


r1  k-i  „  J  '2  „ 

f  u  ^  (1-u)  du 

0 


r(k  +|)  r(|  ■  1) 

r(k  +  §  -  i) 


Thus 

r(i)  r(£  -  1) 

Iq  =  r(|  -  i) 

■f*  r(§  -  1) 

r(2-|-i) 

and 

T  _  1.5.5.  ( 2k- l) 

k  “  (n-l)(n+l)  *.»  (n+2k-3j 


(2.2.37) 
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Substituting  equation  (2.2.37)  in  equation  (2.2.35)  we  get 


23 

n( l-p2)2( 1-r2)  2 


2it(  l-2pr+p2)2 


00 


k=l 


cp(^k)(r) 


_ LLi _ __LL1 

[(n-l)(n+l)  ••• 


iRk-DL 


(n+2k-3) ] 


or 


(2.2.38) 


n  r(|  -  l)(l-p2)2(l-r2)2 
h(r ) - § - - 

2  Vtc  I* (^~)  ( 1-pr)  ( l-2pr+p2)2 


{1  +  0  ( n  " 1 ) ) 


since 


<p(0 


linfii  r  i 

1  -  pr  l  n 


l-2pr+p2 

1  -  pr 


_  U-rg)^ 

1  -  pr 


{1  +  0( n- 1 ) 3  . 
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•  i:  ton 


([23]  in  Chapter  II  Section  4)  has  shown  the  following; 


(i) 


r(|  -  1) 


is  0(n  2) 


(ii) 


h(p) 


n 


{  2jt(l  -  P2) 


1 

2 


(1  +  0(n-1)) 


(iii)  h(r)  is  0(n2)  at  r  =  p  , 

(iv)  Var  (r)  is  0(n  which  suggests  that  the  distribution  of  r 

_  1 

is  concentrated  about  r  =  p  within  a  range  which  is  0(n  2) 
on  either  side  of  p  . 

_  JL 

(v)  Over  the  effective  range  of  r,  (r-p)  can  be  taken  as  0(n  2). 

(vi)  If  r  is  replaced  by  p  in  the  first  neglected  term  in  the 

expansion  of  h(r)  the  term  is  altered  by  an  amount  0(n  . 


(vii) 


11  -  P2)2(l  ~  r2)2  _ 

1  -  Pr 


l-2Pr+p2  i 


2(1-P2) 


[1  +  0(n"'^2)] 


(viii) 


N-l 


h(r)  K 


JLkrfI 


(l-2pr+p2) 


N 

2\2 


[1  +  0(n’5/2)] 


where  k  is  an  adjusted  normalising  constant  and 


N  =  n  -  1  + 


1  -  P‘ 
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Below  we  describe  the  renormalization  of  h(r)  as  detailed  by 
Patton  [25]  Chapter  II  Section  4. 


Consider 


k 


-1 


k-1 

1  (l-r2)  2 

k 

(l-2pr+p2)2 


dr 


k  =  1,2,.. 


and  let 


u  = 


1  +  r 


Then 


Jk  = 


(W 


k-1  k-1 


/V  (1-u) 
J  0 


1  _  ^ 

(l+p)2 


k 

2 


du 


It  can  be  shown  (see  [41 ],  page  293>  example  No,  l)  that 


(2.2.39) 


J.  = 


2k  [r(i^l)]2 


(1+P) 


k  r(k+i) 


k  k+l 
2  ’  2 


;  k+l  ; 


-ie. 


(1+p)' 


where  F(a,b;c;z)  is  the  usual  hypergeometric  function.  It  follows  from 
([10],  page  64,  formula  24)  and  the  hypergoemetric  power  series  that 


k  k+l  4p 

2  ,  2  ;  k+l  ;  (1+p)2 


=  (1+p)  F 


k  .  k  1  o 
-  ,  0;  -  +1;  p2 


=  (1+p) 


Then  equation  (2.2.39)  becomes 


2k  [r(is|i)]3 


r(k+i) 


-  58  - 


Using  the  duplication  formula  for  the  gamma  function,  ([4l],  page  240 ) 
have 

2k  ir(Jal)  ]2 

Jk  =  2k  (*)■*  r(lal)  r(|+i) 

r(*fh 
r(|+  i) 

Thus  for  k  =  1,2,...  we  have 
(2.2.40)  Jk  =  f 


JLklf) 


k- 1 
2\  2 


dr  = 


( l-2pr+p2) 


2\  2 


r(ial) 

r(|  + 1) 


Employing  this  general  result,  we  see  that 


I  = 


P1  ( l-ra) 


J 


-1 


N^l 

2 


( l-2pr+P2) 


dr 


r(^) 

r(|  +  i) 


Then  the  renormalized  density  function  is  given  by 

N-l 

2 

(2.2.41)  h(r)  ^  K(n,p)  ■^1~r  ^ - -  [1  +  0(n“^/2)J 

( l-2pr+p2)2 


we 


■  -j  -if  )  ,  r,  •  ■'  c  -■  i  i  -J-.  i*  -  ■  i 
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where 


(2.2.42) 


K(n,p)  =  7  = 


r(f+i) 
/*  r(S|l) 


and 


N  =  n  -  1  + 


1  -  P‘ 


This  is  Leipnik's  [20]  form  of  h(r)  . 

Daniels  [5]  also  considered  the  case  of  the  Unknown  Mean, 
we  summarize  his  results. 

When  the  mean  is  unknown,  it  is  estimated  by 

5  X.  +  X_  +  .  •  .  +  X  ,+ix 

^  1  2  n-1  *  n 

x  =  - - - - - 

n  -  1 

The  serial  correlation  coefficient,  p  ,  is  estimated  by 

C 

r  =  “ 
o 


where 

C  =  (x1-x)(x2-x)+...+(xn_1-x)(xn-x)  =  c- (n-l)x2 

and  C  =  i(x  -x)2+(x  -x)2+. . ,+(x  ,-x)2+^(x  -x)2  =  c  -(n-l)x2 

o  ^  1  '  v  2  n-1  '  *  n  '  o 

with  c>c0  defined  as  in  equations  (2.2.6)  and  (2.2.7). 


Below 


1|0  - 


where  C 


with  B 


and 


we  have 


The  joint  moment-generating  function  of  C  and  Cq  is 


T  C  +  TC 

M(T  ,T)  =  £(e  °  °  ) 


.  ikefl2 


n 


(2ir)2 


poo  T  C  +TC-ix'C  ^x 
...  /  e  °  °  dx  ...dx 

J  1  n 


-00  -00 


is  given  by  equation  (2.2.3).  Since 


T^+TC-^x'C” lx  =  Tnc^+Tc-|x'c"1x-(n-l)(To+T)x2 


o  o 


=  "2  x'3  X 


(T  +T) 


x'  mm'  x 


n  -  1  - - 


r  2(T  +T) 
i  x'  |  B  + - °’Y  52  s’ 


x 


given  by  equation  (2.2.8), 


X  —  (xt,  •••»  x^) 


m'  =  U,  1,  1,  1,  1,  J)  , 


M(Tq,T)  = 


IkPfi2 


n 


poo  p  oo 

J  "  J 


-fe' 


2(T  +T) 

„  o 

B+  — ™ “  m  in 


x 


(2n) 


2  "CO 


■00 


dx  dx0, . ,dx  , 
12  n 


which  gives  us 


-  hi  - 


M(TQ,T)  =  (l-p2)2  | B  + 


2(T  +T) 
o 

- : —  m  m 

n  -  1  —  — 


_l 

'2 


2(VT) 

On  evaluating  |B  +  - - —  m  m' |  and  making  the  necessary 

substitutions  and  reductions,  we  obtain  the  following  approximation  after 

. , . .  .  .  n- 1 

omitting  a  term  m  z  : 


n- 1 


M(u-rT,T)  ^ 


( l-p)2( l-z2)^( 1-z) ( l-2rz+zs)  2 


(l-p)[ ( 1-pz) [ l-pr-(r-p)z]-u( 1-z2)] (l-2pr+p2-2u) 


n-5 
2  * 


with 


:  jl  + 


( l+p)2( l-£rz+z2) [ ( z-p) ( 1-pz) ( 1-r )+u( 1-z2) ] 
(n-l) ( 1-z2) ( l-2pr+p2-2u) [ 1-pr- (r-p)z-u( 1-z2) ] 


1 


z  +  — 

z 


1  +  P2  -  2T 
_ _ o 

P  +  T 


To  obtain  the  approximation  with  remainder  relatively  0(n  ), 

we  ignore  the  last  factor  and  we  take  the  dominant  term  as  before  in  the 
expansion  of  the  integral  for  h(r).  On  renormalizing,  we  have 


n^ 

2 


h(r)  „  K(i-og)gd,a  n  +  0(n-3/a)] 

( 1-pr) ( l-2pr+p2)  2 


where  K  is  the  normalizing  constant.  In  Leipnik's  form, 


-  42  - 


(2.2.43)  h(r)  ^  ~ 


2  r  (2*2) 


I'1 


r  (|)  [n( i-p)+( i+p) ] 


(l-r)(1~r3)  U  +  0(n'5/2)]  , 


( l-2pr+p^) 


N- 1 

2\  2 


with 


N  -  n  - 1  + 


3.  Non-circular  Markov  Process:  Patton's  Case. 

Consider  a  non-circular  Markov  Process  described  in  (2.1,1)  as 
taken  by  Daniels,  Take  the  sample  estimate  of  p  to  be 


r 


where 


c 


1 

2 


xi  +  xixa  + 


+  X 


n-1 


x2 

n 


(2.3.1) 


22  2 

c  =  xf  +  x^  +  . . .  +  x^ 
o  1  2  n 


Following  Daniels1  Method,  Patton  [25]  obtained  the  following: 


The  joint  moment-generating  function  of  c  and  cq  is 


T  c  +Tc 

M(To>T)  =  £(e  °  °  ) 


( 1-Pa)~2 

(80* 


-00 


T  c 
o  o 


+Tc-J  x'  C_1~ 


dx. . . . dx 

1  i 


-00 
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where  C  is  given  by  the  equation  (2.2.3). 


Since 


ToVTc  -is’o'h"  -i  [(1-2To-T)x2  +  (i+p2.2Io)(x|+...+x*_1) 


where 


+  (1-2Tq-T)x2  -  2(p+T)(x1x2+,..+xn_1xn)] 


^  x*  D  x  , 


(2.3.2)  D 


1-2T  -T  -(p+T) 


(P+T)  l+p2-2To  -(p+T) 


(p+T)  1+P2-2Tq  -(p+T) 


(0) 


(0) 


(p+T)  l+p2-2T 


-(p+T) 


and 


x 


(xx,  • • • »  XR) 


it  is  seen  that 


(2.3.3)  M(T  ,T) 


=  Ik £ 

(ait ): 


-ns)2  p  J 


\00  l 


x'D  x 


dx. dx  . , . dx 
12  n 


-CO 


-(0+T) 
1-2T  -I 


-oo 


Diagonalizing  D 


and  integrating,  we  obtain 


(2.3.U)  m(Tq,t)  =  (l-p2)2  |d|'2 

_1 

where  |  D  j  2  is  the  Jacobian  of  the  transformation. 

On  evaluating  jD|  and  making  necessary  substitutions  and 
simplifications  it  is  found  that 


,n-2 


(2.3.5)  |  D  |  ^  .(l-gpr+p^-gu),  [  ( 1-pz)  ( l-2pr+pz+p-z)-2u(  1-z)  ]; 

( 1-z2) ( l-2rz+z2)n 


and  therefore 


M(u-rT,T) 


( 1-P2)2( l-z2)2( l-2rz+z2) 


H-l 

,2 


( l-2pr+p2-2u)  [ ( 1-pz) ( l-2pr+pz+p-z)-2u( 1-z) ] 


Then  using  (1.2)  we  obtain 


(2.3.6)  h(r) 


n( 1-P2)2 


2jti(  l-2pr+p2) 


2  ■ 1 


f  cp(z)  (l-2rz+z2) 


—  -  2 

2 


dz 


where 


Uizfl 


3/S 


(2.3.7)  <p(z)  -  ( 1-pz) ( !-2pr+pz+p 


(1  + 


n 


2(  1-z) ( l-2pr+p2) 

( 1-pz) ( l-2pr+pz+p-z) 


-  k 


To  obtain  an  approximation  with  remainder  relatively  0(n  J'  ~) 
the  last  factor  was  ignored  and  the  dominant  term  was  taken  in  the  expansion 
of  the  integral  for  h(r). 


-  1*5  - 


On  renormalizing  h(r)  was  obtained  as 

N+l 

(2.3.8)  h(r)~K(n,o)  - - ;j  [l+0(n'3/2)] 

(l-r)(l-2pr+p2)^ 

where 

2  r(|  +  2) 

(2.3.9)  K(n,p)  =  - - - 

sT*  r(^)  [N(  1+p)  +  2] 

and 


2 

(2.3.10)  N  =  n  -  1  +  -f-p2 


which  is  in  Leipnik's  Form. 


If  the  mean  is  unknown  we  take  as  the  estimate  of  the  unknovm 


mean 


(2.3.11) 


x  = 


X,  +  ...  +  x 

1_ n 

n 


and  the  estimate  r  of  the  serial  correlation  coefficient,  p  ,  to  be 


(2.3.12) 


r 


C 

C 


where 


(2.3.13) 


r 

C  =  |(x-x)2  +  (xx-x)(x2-x)  +  ...+  (xn_rx)(xn-x)  +  i(xn-x)2 


4 


C 

o 


c  -  n  x2  , 

<v*)2  +  •••  +  <v*)S 


=  c  -  n  x2 
o 


)-i6  - 


with  c  and  as  defined  in(2,3.1).  Using  Daniels'  method  it  is 

that 


(2.3. 1U) 


N 

h(r)  ^  .K  -  [1  +  0(n"5/2)] 

N-l 

( l-2pr+p2)  2 


where  K  is  an  adjusted  normalizing  constant  and 


N  =  n  -  1  +  Y  -  p'2 


In  this  case,  the  renormalized  density  function,  h(r),  in  Leipnik1 
is  obtained  as 


(2.3.15) 


h(r)  ^  K(n,p) 


Ikrfl 


( l-2pr+p2) 


N-l 
2\  2 


[1  +  0(n"3/2)] 


where 


2  r(5t2) 

(2.3.16)  K(n,p)  =  — - * - 

r(|+l)  [(l-p2)(N-l)+lt] 


and 


found 


s  form 


b7 


CHAPTER  III 

THE  APPROXIMATE  DISTRIBUTION  OF  ESTIMATES  OF  THE 
SERIAL  CORRELATION  COEFFICIENT  FOR  A  PRESTATIONARY 
LINEAR  MARKOV  PROCESS  WITH  KNOWN  MEAN 


In  the  process  defined  by  (2.1.2)  following  Daniels  [5]  we  take 


the  sample  estimate  of  p  to  be 


c 

c 

o 


where 


(3.1.1) 


c 

=  xlx2 

+ 

fT\ 

M 

CM 

X 

4“  •  •  • 

r-H 

I 

+ 

9 

c 

0 

„  JL  v2 
-  2  x! 

+ 

2 

X2  + 

xf;  4* 

j 

...  +  x2  . 

n- 1 

+  i  x2  , 
^  n  " 

and  employ  a  method  similar  to  that  used  by  Daniels  for  obtaining  the 
approximate  distribution  of  r. 

In  matrix  notation  (2.1.2)  takes  the  form 


x  =  R  e 


where 


«■ 

■■■  — 

X1 

1 

0 

0 

•  •  • 

0 

0 

el 

X2 

p 

1 

0 

•  •  • 

0 

0 

e2 

• 

• 

,  R  = 

P2 

P 

1 

•  •  • 

0 

0 

and  e  = 

• 

• 

• 

• 

• 

• 

#  •  • 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

• 

X 

n 

*n-l 

P 

n-2 

P 

d"'3 

•  •  • 

•  •  • 

• 

p 

• 

1 

e 

n 

Hence 


e 
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We  can  also  write  (2.1.2)  in  the  following  form 


e 


1 


=  x 


1 


e 


2 


e 

n 


x 

n 


pv 


1 


which  gives  us 


— 1 

J— 4 

_ 1 

X, 

-0  1  (0) 

1 

X2 

-p  1 

X, 

3 

• 

(0)  -p  1 

-p  1 

• 

• 

X 

n 

Hence 


and 


' 
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(R"1)'  - 


1  ~P 

1  “P 

1  ~P 


(0) 


(0) 


1  -p 
1 


Thus 


C"1  =  (R-1)'  R_1 


(3-1.2) 


1+P2  -p 

-p  1+p2  -p  (0) 

-p  1+p2  -p 


(0) 


“p 


1+p2  -p 
-p  1 


Hence 


x'  C  1  x  =  x*  (R  ■*■)'  (R  x 


=  [ x^ ( 1+p 2)  -  px2,  -  px^  +  ( 1+p2) Xp  -  px^ 


px  +  ( 1+p2) x  -  px  ,  -  px  +  x  ] 

n-2  n-1  n  n-1  n 


x. 


x. 


n 
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=  x2(l+p2)  -  p  x1x2 

+  x|(l+p2)  -  0  x^  -  p  x2x^ 

+  x2(  1+p2)  -  p  xgx  -  p  x  x^ 

+  x2(  1+p2)  -  p  x^  -  p  x^x^ 


+  Xn-2^1+°2^ 
+  x2_^(  1+p2) 


P  x  _x 
n-3  n-2 


P  X  -X  . 
n-2  n-1 


P  x  0x  . 
n-2  n-1 


P  X  -X 
n-1  n 


-  P  x  -X 
n-1  n 


n-1 


=  ( i+P3)  £ 


-  2p 


i=l 


n-1 

I 

i=l 


x.x.  ,  +  xc 
i  i+1  n 


The  joint  distribution  of  x. ,  x_,  ...»  x  is  given  by 

12  n 

n 

*>  p  _  X  _  n 

dF(x1,x2,...,xn)  =  (2*)  f$|  2  exp(-Jx'  0.  *)  dx]dx2...dxn  . 

Now 

|c-1|  =  Kr"1)'  r"1! 

-  I*'1!2 

=  i  . 
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Thus 


n 


dF(x1,x2> . . .  ,xn)  =  (2«)  c  exp(  x'  C  1  x)  dx^^dx^.^dx^ 


Then  the  joint  moment-generating  function  of  c  and  cq  is 


M(To>t)  =  £[exp(Toco+Tc)] 


n 


■  r-s 


>oo  T  c  +Tc  — \  x' C  ^  x 

e  dx,dx  . ,.dx 

12  n 


-CO  -00 


-1 


where  C  is  given  by  (3.1.2),  Now, 


T  c  +Tc-  \  x'£-1  x  =  -  2 


o  o 


n-1  n- 1 

{(i^)7  x^  ~  2p  y 

ill  i=l 


2p  )  x.x.  .  +  x‘ 
i  l+l  n 

i=l 


-  2T 0(ix2+x|+...+x2_1+ix2) 

-  2T ( Xlx2+x2x5+ ...  +xn_ 1xn )  j* 

i  |(l+P2-To)x2+(l+p2.2To)(x|+...+x2_l)+(l-To)x2 

-  2(p+T)(x1x2+x2x5+.. .+xn_1xn)  j- 


lx'  B  x 


where 


' 
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i+p2-t 

o 

-(P+T) 


B  = 

(5.1.3) 


"(P+T) 

1+P2-2T 

o 

"(P+T) 


(0) 


"(P+T) 

1+P2-2To  -(p+T) 


-(P+T) 


(0) 


1+p2- 2T  -(P+T) 
o 

-(p+T)  1-T 

o 


Then 

_  n 

~Z  poo  p  00  _i  , 

M(T  ,T)  =  (2k)  /  ...  /  e"2---  dx....dx 

o  J  J  In 

-00  -00 


As  mentioned  in  Chapter  II,  Section  2  ,  there  is  an  orthogonal 

transformation  of  that  will  diagonalize  B.  Diagonalizing 

B  and  carrying  out  the  integration,  we  obtain 

(3.1.4)  m(to,t)  =  |B|"S 

where  |b|  2  is  the  Jacobian  of  the  transformation. 

| B |  we  proceed  as  follows;  (p+T)  is  factored 
and  the  following  substitutions  are  made: 

1 

"P  +  T 

1  +  P2  -  2T 

o 

^  P  +  T 

so  that  | B |  may  be  written  as 


To  evaluate 


from  each  row  of  |Bj 


(3-1.5) 


■i 


»■ 
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(3.1.6) 


I —I  =  (P  +  T)n  Dn 


where  D  is  the  determinant 
n 


D  = 
n 


b+aT 


b  -1 


-1 


(0) 


b 

-1 


-1 

b 


-1 


(0) 


-1  b  -1 


-1  a( 1-T  ) 
o' 


n 


b  -1 


•1  b  -1 


(0) 


(0) 


•1  b  -1 


■1  a(l-T  ) 
'  o ' 


n 


aT 


-1 

b 


(0) 


-1 


■1  b  -1 


(0) 


b  -1 


•1  a( 1-T  ) 
o 


n 
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Or 


(3.1.7) 


D  =  A  +  a  T  A., 
n  n  o  n-1 


where 


(3.1.8)  A  = 


n 


b  -1 


-1  b  -1 


(0) 


(0) 


-1  b  -1 


-1  b  -1 


-1  b  -1 


-1  a( 1-T  ) 
'  o' 


n 


and  A-n_^  is  obtained  from  A^  by  removing  the  first  row  and  first 
column. 


Therefore 


(3.1.9) 


i 


‘1  -  a<X  -  To>  > 


A  =  ab(l  -  T  )  -  1  , 
2  o' 


and  in  general. 


A.  =  b  A.  -  A. 

3  J-l  J-2 


or  equivalently 


(3.1.10) 


A.  -bA.  +  A .  =0 
3+2  J+l  J 
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Taking  E  to  be  the  forward  difference  operator  of  the  calculus 
of  finite  differences  defined  by  equation  E(A^)  =  ^j+i  *  we  ma^  wrlte 
equation  (3.1.10)  as 


(E2  -  bE  +  1)A  =  0 


Solving  this  difference  equation,  we  find  that  the  roots  of  the 
auxiliary  equation 

z2  -  bz  +  1  =  0 


are 

b  +  \fb2  -  4 
2  *  - 2 - 

and 

1  _  b  -  s/b2  -  4 

z  2 


where 


(3.1.11)  z  +  7  =  b 


Hence  the  general  solution  of  equation  (3.I.I0)  is 
(3.1.12)  A.  =  klZJ  +  k,,z'j  . 

Applying  conditions  (3.I.9)  to  equation  (3. 1.12),  we  have 

Al  =  a(l-TQ)  =  kxz  +  k2z  1  , 

A2  *  -  1  =  V2  +  V'2  • 

From  these  equations  and  equation  (3.I.H)  the  values  of  k^  and  k^  are 


found  to  be 
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1  -  a(l-T  )z 
o 


^a(1 T0)  "  z^z 


1  -  z‘ 


1  -  z‘ 


Substituting  these  values  in  equation  (3. 1.12),  and  taking  j  =  n  we  have 


(3.1.13) 


or 


(3.1.14) 


where 


f l-a( 1-T  )z^ 
a  J  o'  I  n 

'  {a( 1-Tq)-z)z  - 

A  -  1  r  z  + 

n  ^  1  -  z2  [ 

2 

L  1  -  z 

n  p 

-n 


A  =  •—  5 
n  1-z^ 


(1-z2) 


—  -  a(  1-T  )  —-5  ( 
n  '  o'  1-z^  ' 


z  ,  n 

z  -  z 


n)  . 


a  = 


b  = 


P  +  T 

1  +  P2  -  2T 
_ c 

P  +  T 


1 

z  +  - 
z 


Substituting  equations  (2.2.18),  (2.2.20)  for  u  and  (p  +  T) 
respectively,  we  see  that 


a  = 


1  -  2rz  +  z2 _ 

z( 1  +  p2  -  2pr  -  2u) 


1  +  P2  -  2T  =  (^t  1)  <  V  0*  -  SB*  -  2Ml 

o  v  '  1  -  2rz  +  z^ 


From  equation  (2.2.22) 


(3.1.15)  1  -  To  =  1  -  j 


( l+z2)(pr+u)-rz( 1+p2) 

l-2rz+z2 


l-2rz+zg-( 1+z2) (pr+u)+rz( 1+p2) 

l-2rz+z5 
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Hence 


a(l-T  )  = 
'  o ' 


l-2rz+z2-(  l+z2)(pr+u)+rz(  1+p2) 

z( l+p2-2pr-2u) 


( l+z2)( l-u-pr)-rz( 1-p2) 

z( l+p2-2pr-2u) 


Then  equation  (3. 1.11+)  becomes 


A 

n 


z2n-z2  ( z^H-l)  f  ( 1+z2) ( l-u-pr)-rz( 1-p5) 
z"(l-z2)’  (1-z2)  z"  l  ( l+p2-2or-2u) 


or 


(3.1.16)  A 

n 


1 

zn(l-z2) 


’  (l-z2n)  [■(  1+Z2)(  1-u-pr )-rz(  1-p2)) 

l+p2-2pr-2u 


?(1' 


2n-2\ 
■z  ) 


Again  using  equations  (3. 1.14)  and  (3.I.I5)  we  have 


(3.1. !7)  aTo  An_1 


f ( 1+z2) (pr+u)-rz( 1+p2) ]  "( 1-Z2n  2) f ( 1+z2) ( 1-u-pr )-rz( 1-p2) 

z( l+p2-2pr-2u)zn  ^(l-z2)  L  l+p2-2pr-2u 


2/ ,  2n-l+v 
-  z^l-z  ) 


Hence  using  equations  (3.1. 5),  (3.1.6),  (3. 1.14)  and  (3.I.I5)  we  have 


1*1  ( l+p2-2pr-2u)n  1 _ 1 _ 

"  =  ( l-2rz+z2)n  l+04-20r-2u 


( l-z?n){ ( 1+z2)  ( 1-u-pr )-rz( l-o2) } 


2n-4 


z2(l-z^n  ^){ ( 1+z2) (pr+u)-rz( 1+p2) }-  z2(l-z2n  2) ( l+p2-2pr-2u) 


( l-z2n  2)f ( 1+z2) ( 1-u-pr )-rz( 1-p2) If ( 1+z2) (pr+u)-rz( 1+p2) } 

1+P2  -  2pr  -  2u 


x 
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2n  2n-2 
z 


2n-k 


We  omit  the  terms  in  z~  ,  z  and  z  within  the 

bracket  which  gives  rise  to  an  error  that  is  exponentially  small.  Hence 
ignoring  these  terms  we  get 


n- 1 


B 


(l+pg-2pr-2u) 

( l-2rz+zg)n( 1-z2) 


( l+z2)( 1-u-pr )-rz( 1-p2) 


z2{ ( 1+z2) (pr+u)-rz( 1+p2) } 


-  z2  ( l+pg-2pr-2u) 

{ ( 1+  z  2 )  ( l-u-pr)-rz( 1-p2) 1 f ( l+z2)(pr+u)-rz( l+P2) 1 

+  1  +  pa  -  2pr  -  2u 


( l+ps-2pr-2u)n~2 

( l-2rz+z2)n( 1-z2) 


{ ( 1+z2) ( 1-u-pr ) -rz( 1-p2) ){ l+p2-2pr-2u] 


-  z2{pr+u+prz2+uz2“rz-p2rz+l+p2-2pr-2u] { l+p2-2pr-2u] 

+  { ( 1+z2) ( 1-u-pr )-rz( 1-p2) }[ ( l+z2)(pr+u)-rz( 1+p2) ) 


( l+p2-2pr-2u)n  2 

( l-2rz+z2)n( 1-z2) 


1+z2) ( 1-u-pr )-rz( 1-p2) ) 

* { l+p2-2pr-2u+pr+u+prz2+uz2-rz-p2rz ) 

-  z2{  l+p2-2pr-2u]  { l+p2-pr-u-rz(  l+p2)+z2(pr+u)) 


( l+p2-2pr-2u)n  P( l+p2-pr-u-rz( l+p2)+z2(pr+u) } 

(l-2rz+z2)n  (1-z2) 


l-u-pr-rz( l-p2)+z2(u+pr-p 


X 
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or 

(3.1.18)  | B |  ^  (l+P2-2pr72u)n’2  (l+P2-p.f-u-rz(_l+P2>zf(£r+u.)] 

(l-2rz+z2)n  (1-z2) 


»  { l-u-pr-rz( l-p2)+z2(u+pr-p2) ) 


Thus  an  approximation  of  M(Tq,T),  as  given  by  equation  (3. l.U) , 


is 


(3.1.19)  M(u-rT,T)  ^ 


(1-z2)2  ( l-2rz+z2) 


n-2 

2 


( l+p2-2pr-2u)  (l+p2-pr-u-rz(l+p2)+z2(pr+u) )2 


{ l-u-pr-rz( l-p2)+z2(u+pr-p2) )2 


Combining  equations  (2.2.25)  and  (3.I.I9)  we  get 


M(u-rT,T)  - tV2»^2 


S-2 


( l+p2-2pr-2u)  2  { l+p2-pr-u-rz( l+p2)+z2(pr+u) )2 


* 


1 

- T 

( l-u-pr-rz( l-p2)+z2(u+pr-p2))2 


Differentiating  partially  with  respect  to  u,  we  get 


-  6o  - 


a 


M( 


n  g 

u-rT,T)  ~  (1-z2)5/2  (l-2rz+z2)2  ( l+p2-2pr-2u) 

^  _(n~  h.\ 

+  i-  { l+p2-pr-u-rz(l+p2)+z2(pr+u) }  ^  Cp  2  <p”2 

uu  J  d 


(inks 

2  ’  -i  -1 


+  ^  { l-u-pr-rz(  l-p2)+z2(u+pr-p2) )  2  op^  2  cp^2 


where 


cpj^  =  1  +  p2  -  pr  -  u  -  rz(  1+p2)  +  z2(pr+u)  , 


cp2  =  1  -  u  -  pr  -  rz(  1-p2)  +  z2(u+pr-p2)  , 


©  =  1+p2  -  2pr  -  2u  . 

3 


Thus 


a 


M(u-rT,T)  |f 


(1-z2)5/2  (l-2rz+z2)2  ij>j  2  (-2)  oj2  <pj 


3,  '  1 

" 2 


-  i  <P^5//2  (-1+22)  l>= 


-(^)  1 

2  «P„2 


-  2  ("1+z2)  IP5  2  ‘Pi2 


5-2  1  .1  r 


(l-z2)3//2(l-2rz+z2)2  q>j2  (p22  <P5  2  n-U  +  5(l-z2)cpj(<p11+q>21) 


-  61  - 


Hence 


a 


M(u-rT.T)  |f 


u=o 


n(l-z2)3/2  ( l-2rz+z2| 


?  -  2 


-(— ) 

( l+p2-2pr )  2 


{  l+p2-pr-rz(  l+p2)+prz2}  +  2{  l-pr-rz( l-p2)+(pr-p2)z2)+2 


A 


u 

n 


+  JL  ( l-z2)( l+p2-2pr)(2-2pr+p5-2rz+2pr-p2z2) _ 

2x1  { l+p2-pr-rz( l+p2)+prz2){ l-pr-rz(  l-p2)+(pr-p2)z2) 


Thus  by  equation  (1.2),  we  have 


(3.1.20) 


h(r) 


n 


2«i(  l+p2-2pr ) 


n-2 

2 


n 

J'  <p(z)  (l-2rz+z2)2 


dz 


where 


(3.1.21)  cp(z) 


(1  -  Z2) 


5/2 


{ l+p2-pr-rz( l+p2)+prz2)2( l-pr-rz(l-p2)+(pr-p2)z2}2 


1  r(l-z2)(l+p2-2pr)(l-pr+|  -rz+pr-|  z2) 

1  +  1  iiy -  - — 

n  l 


{ l+p2-pr-rz(l+p2)+prz2}{  l-pr-rz(  l-p2)4-(pr-p2)z2} 
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As  shown  in  Chapter  II,  Section  2,  the  transformation 


z  +  ± 
z 


1  +  P2  4-  2rT 

P  +  T 


maps  the  T-plane  cut  exterior  to  the  interval 


-( 1+p)2  ( 1-p)2 

2(l+r)  ’  2(l-r) 


on 


the  real  axis  onto  the  interior  of  the  unit  circle,  |z|  =  1,  in  the 


z-plane.  Also  the  path  of  integration  crosses  the  real  axis  at  z  =  p 

-  id  id 

and  terminates  on  the  boundary  of  the  unit  circle  at  e  and  e  , 


where  r  =  cos 9. 


The  only  possible 
of  h(r)  occur  if 

either  cp^(  z)  = 

or  q>2(z)  = 

Discussion  of  cp^(z) 


singularities  of  the  integrand  of  the  integral 

1  +  p2  -  pr  -  r(l+p2)z  +  prz2  =  0  , 

1  -  pr  -  r( l-p2)z  +  P(r-p)z2  =  0  . 


Nov? 


cp1(r)  =  1  +  p2  -  pr  -  r2(  1+p2)  +  pr^ 

=  ( l+p2)( l-r2)-pr( 1-r2) 

=  ( 1-r2) ( l+p2-pr ) 

>  0  7  since  1  -  r2  >  0  and  1  -  pr  >  0  implies 

1  +  p2  -  pr  >  0 

Cp1  ( p)  =  1  +  p2  -  pr  -  r(l+p2)p  +  prp2 
=  1  +  P2  -  2pr 
=  (1-pr)2  +  P2(l-r2) 
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>  0  since  (l-pr)2  >0,  p2  >  0,  1  -  r2  >  0  for 
| r  |  <  1  .  Now  cp|(z)  =  -r( 1+p2)  +  2prz  =  0  if  z  =  -  . 

|  =  —  >  1  .  Hence  (pj(z)  =  0  for  |z|  >1  . 


But 


2p 


Thus  cp^(z)  does  not  vanish  between  p  and  r  as  |p|  <  1. 


r  <  1. 


Since  Cp^(r)  >  0,  cp^(p)  >  0  and  Cpj(z)  /  0  for  |z|  <  1, 
we  conclude  that  there  is  no  real  singularity  between  p  and  r. 


Now  if 


cp^(z)  =  prz2  -  r(l+p2)z  +  1  +  p2  -  pr  =  0  , 


z  = 


r(l+p2)  ±  \fr2(  1+p2)2  -  lipr ( l+p2-pr) 
2pr 


1+Pg  -J r2(l+p2)2  -  4pr(l+p2-prT 
2p  1  2pr 


The  roots  will  be  imaginary  if  r2(l+p2)2  -  Upr(l+p2-pr)  <  0  .  In  this 

>  1  .  Hence  even  if  there  is  a  complex 


case  Re  z  = 


1 1+p2l 

1+ 

P 

2 

!  2p  ! 

"  2 

P 

root  it  lies  outside  the  unit  circle.  We,  therefore,  conclude  that  cp^(z) 
does  not  vanish  for  any  z  in  the  closed  interval  between  p  and  r. 
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Discussion  of  cpg(z) 

Now  ({^(r)  =  1  -  pr  -  r(l-p2)r  +  p(r-p)r2 

=  1  -  pr  -  r^  +  r  p 

=  ( 1-pr) ( 1-r2) 

>  0  ,  since  1  -  pr  >  0  and  1  -  r2  >  0  . 

<P?(p)  =  1-pr  -  r  ( l-p2)p  +  p(r-p)p2 

=  1  -  pr  -  pr  +  p^r  +  p^r  -  p 
=  1  -  p  -  2pr  +  2p^r 
=  ( 1-p2)  ( l+p2-2pr) 


>  0,  since  1-p2 

>  0  and  1  +  p2  - 

2pr  >  0 

for  | r |  <  1, 

|P|  <  1  • 

Now 

*2{z)  = 

-r( 1-p2)  +  2p(r-p)z  =0  if 

z  _  r.(i.“pf) 

2p(r-p) 

and 

<p2(z)  = 

2p(r  -  p)  . 

If 

=  0  and 

Cp” ( z )  <  0,  then  there  is  a 

maximum  between  p 

and  r 

Hence  there  would  not  be  a  singularity  between 

p  and  r  due  to 

<P2(z) 

as 

<P2(r) 

>  0  and 

92(p)  >  0  . 

There  might  however  be  a  singularity  if  there  were  a  minimum  for 


z  between  p  and  r,  which  is  so 
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if 

cp^(z)  >  0 

or  if 

p(r-p)  >  0 

or  if  either 

P  >  0,  r-p>0  or  p  <  0,  r-p<0 

Hence  there  could  be  a  minimum  if  r>p>0  or  r<p<0. 


Case  I  r  >  p  >  0 

cpg(p)  = 

2p(r-p)p  -  r ( 1-p2) 

2p2r  -  2p^  -  r  +  rp2 

-  r  -  2p^  +  rp2  +  P . 2pr 

< 

-  r  -  2p^  +  rp2  +  p(l+p2)  [Because  l+pa>2p>2pr] 

2  3 

-  r  +  rp^  +  p  -  p^ 

(p  -  r ) ( 1  -  p2) 

< 

0  since  p  <  r  and  1  -  p2  >  0 

And  cp^(r)  = 

2p(r  -  p)r  -  r ( 1  -  p2) 

2pr2  -  2p2r  -  r  +  rp2 

-  r( 1  +  p2  -  2pr ) 

< 

0  since  r  >  0  ,  1  +  p2  -  2pr  >  0 

Hence  ( p^(z)  is  ne 

gative  for  z  =  p  and  z  =  r  . 

is  negative  for  z  =  p  and  z  =  r 
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Case  II  r  <  p  <  0 


Let  r '  =  -r  and  p '  =  -p 


Now 


<Pg(p)  = 


2  3 

r  +  3rp  -  2p 


=  r'  -  3r'p'2  +  2p ' 


=  r'  -  r'p'2  +  p ' ( -2r 8  p ' )  +  2p 


,3 


>  r'  -  r'p'2  +  p'{  -(l+p'2)}  +  2p'^ 
=  r'(l-p'2)  -  p '  -  p8^  +  2p ' ^ 

=  (r 8 -p 8 )  ( 1-p 8 2) 


£  since  1+p*2  >  2p8r8, 
-(l+p'2)  <  -  2p 8  r 8 


>  0 


p 

since  1-p*  >0  and  r'  -  p'  >  0  . 


Again  Cp^(r)  =  2pr  -  r  -  rp' 


=  -2p'r’2  +  r8  +  r'p'2 


=  r8  (1  +  p'2  -  2p 8  r 8  ) 


>  0 


since  r8  >  0  ,  1+P1  -  2p8r8  >  0 


,2 


Hence  cp,!>(z)  ‘>  ^  ^or  z  =  p  anc^  z  =  r  . 


Cases  I  and  II  show  that  Cp^(z)  does  not  vanish  between  z  =  p 
and  z  =  r.  Hence  there  cannot  be  a  minimum  between  z  =  p  and  r.  We, 
therefore,  conclude  that  there  is  no  real  singularity  in  the  closed  inter¬ 
val  between  p  and  r  due  to  cp^ ( z )  • 
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Nov/  if 

qp2(z)  =  p(r-p)z2  -  r(l-p2)z  +  1  -  pr  =  0  , 

z  _  +  r(l-p2)  ±\/r2(l-p2)2  -  4p(r-p)(l-pr) 

2p(r  -  p) 

The  roots  will  be  imaginary  if  r2(l-p2)2  -  4p(r-p)( 1-pr)  <  0 
which  is  possible  only  if  p(r-p)  >  0. 

If  there  are  imaginary  roots  let  them  be  x  +  iy.  Obviously 

_  sfb p(r-p)(l-pr)-ra(  1-p^V2 
’  y  "  2p(r-p) 

2  _  ..2  v2  _  Mr-p)(l-pr)  _  1  -  pr 
hp  (r-p)2  p(r-p) 


Hence 

Ms>  i 

if 

Y°\  > 1 

P(r-p) 

or 

if 

1  -  pr  >  p(r-p) 

or 

if 

1  -  pr  -  pr  +  p2 

or 

if 

1  -  2pr  +  p2  >  0 

which  we  know  is  true. 

Hence  the  complex  roots,  if  any,  lie  outside  the  unit  circle. 
We  therefore  conclude  that  there  is  no  singularity  between  p  and  r 
due  to  Cp2(z)  . 


x  = 


rjl-P  2I 

,r-P) 


2p(: 


Therefore 
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Because  there  are  no  singularities  in  the  closed  interval  joining 

p  and  r,  we  may  deform  the  path  of  integration  in  the  z-plane  to  be  the 

•  i  9  i9 

straight  line  joining  e  to  e  and  crossing  the  real  axis  at  z  =  r. 

As  in  Chapter  II,  section  2,  we  may  use  the  substitution  (2.2.29) 
for  z  on  the  path  of  integration  so  that  equation  (3.1.20)  becomes 


h(r)  ^ 


a- 2 

/ 1  2\ 2  2 
n( l-r^J 

n 

2n  (l+P2-2pr)2 


1 


n 

cp(z)  (1-cj2)2 


Following  the  method  of  Daniels  [5  ]  discussed  in  Chapter  II,  we 
expand  cp(z)  as  a  power  series  in  (z-r)  and  integrate  with  respect  to 
u)  .  By  using  equation  (2.2.33) >  it  is  seen  that 


cp( z)  -  cp(r)  +  iu)(l-r2)2  cp'(r)  + 


.k  k,_  2v 
l  a)  (1-r  ) 


kl 


q)(l^(r) 


and  hence, 


(3.1.21)  h(r )  ^  - p-^~r 


2\  2 


—  -  1 


2it  (1-2P-1-P2)' 


00 


k=o 


k 

yfc(r)  J1 


a  , 

2~‘ 


for  the  case  that  k  is  an  odd  integer,  the  integrand  is  an  odd  function, 
and  the  integral  vanishes.  Therefore,  we  may  write  equation  (3.1.21)  as 


(3.1.22) 


00 


r2 


(l-2pr+p2)‘ 


k=o 


-1 


dw  . 


dw  » 


*  . 


- 
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Since 


I  =  r  a)  2  (1  -  a,2  2  da) 
k  J 

-1 


is  the  same  as  equation  (2.2.36),  the  value  of  1^  is  given  by  (2.2.37)* 
Substituting  this  value  in  the  expression  for  h(r),  we  have 


, /  \  n  ( 1-r2) 
h(r)  •'-& - 

( l-2pr+p2) 


l-1 


00 


j  _frx  .  V  ( -  l)k(  l-ra)k  (2k)  [1...?.?, _ ...  (2k-l)i 

V11'  +  /_,  (2k)T  <P  (r)  (n-l)(n+l)...(n+2k-3) 

k=l 


-2_  (1-r8)  2 _  r(2  -1)  „(r)  f  J  +  V 

&/*  1  Tt2^1)  1  ^  2yt)k,'(n-l)(n+l)...(n+2k-3) 


( l-2pr+p2) 


where 


(3.1.23)  <p(r)  = 


(1-r2)2 


( 1-pr )2( l+p2-pr)2 


1  + 


n 


( l+p2-2pr  )( l-pr+“) 
( l+p2-pr)( 1-pr) 


k 


Hence 


Mr) 


n 

2s/t 


(1-r2) 


r(5-i) 


TH 


( l+p2-2pr) 


n-2  p/ti-lv 
2  ^  2  ' 


<P(r) 


i  _  ilnil  SUhl  ,  ofn"2) 

1  2(n-l)  <p(r)  +  °(n  > 


Let 


V'(r) 


( 1-r2)' 


X  i 

( l+p2-pr)2( 1-pr)2 


(3.1.2k) 


- 
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(3.1.25) 

and 


(3.1.26) 


Then 


since 


Or 

(3.1.27) 
Let 

(3.1.28) 
Now 


Hence 


( l+p2-2pr)( 1-pr  +  tt) 

“  ( l+p2-pr)(  1-pr)  14 


z(r)  =  cp(r) 


!  .  Iklfl  2^-1  ,  -2, 

1  at^TT  <p(r)  +  0(n  ' 


z(r)  =^(r) 


1^1- 


-2, 


2n 


<p"(r)  +  0(n  *) 


1 

n-1 


— V  -  Ad 

n(l  -  — ) 

'  n' 


U-l  ln  1  1  V 

n'  nx  n  n^  ' 


=  “  +  0(n“2)  . 

n  '  ' 


z(r)  =  ^(r) 


l  +  q>,’Cr.).)/ 2^.(r)  +  o(n"2) 

n  v  ' 


t(r\  _  /rN  .  l.l..Tr — y.lK.O. 

_  g(r;  2^(r) 


t(r)  =  t(p)  +  t'(p)  (r-p)  +  ^7  t”(p)(r-p)2  +  ... 


=  t(p)  +  t'(p)  0(n  2)  .  [Using  (III. 7)] 


z(r)  =  ^(r) 

=  ^(r) 

=  ctf/(r) 


1  +  C(p)  +  +  0(n-2)' 
n  '  ' 


1  +  +  0(n"3/2) 

n  ' 

1  +  0(n"5/2) 


where  c 


=  1  + 

n 
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Hence 


(3.1.29) 


h(r) 


K 


£  -  1 

(l-r2)2 


-  -  1 

(l+ps-2pr)2  ( l+p2-pr)2( 1-pr)2 


1  +  0(n"5/2) 


where  K  is  an  adjusted  normalizing  constant. 


Let 


(3*1. 30)  f(r)  =  [ ( l+p2-p  r)  ( 1-pr) ]" 2 


Then 


in  f(r)  =  -  \  [in  (l+p2-pr)  +  in(l-pr)]  , 


and 


iilil  =  .if - -  .  -i-p—  1 

f  (r)  2  f  1  +  p2  -  pr  1  -  pr  J 


P  f  l-pr+l+p2-pr 

2  f  fl+p^-pr j ( 1-prJ 


1  p(2-2pr+p2) 

2  ( 1-pr+p2) ( 1-pr) 


(3.1.31) 


f  i  / r\  _  1  p(  2-2pr+p2) 

1  \  ^  *  ~  2  z  /p 

{(l+p2-pr)(l-pr))^ 


Taking  logarithms,  we  get 


in  f'(r)  =  in  |  +  in(2-2pr+p2)  -  ^  in(l+p2-pr)  -  |  in(l-pr) 
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Differentiating  with  respect  to  r,  we  get 


(3.1.32) 


f"(r)  =  -2p  3  P  3  P 

f 1 (r)  2-2pr+p2  +  2  l+p2-or  +  2  1-pr 


From  (3.I.3O),  (3.I.3I)  and  (3.1. 32),  we  get 


f(P)  =  (1-P2)"2 


f'(p) 


p(2-P2) 

2(l-02)5/8 


and 


f"(o)  =  D^-°SX 

2(1-P2)3/2 


~2p  3  3  P 

2-p^  +  2  P  +  2  1-P2 


P2(2-P2) 

"  Ml-P2)5/2 


•H+hp2+6-9P2+3pS-6-3P2 


( 2-p2) ( 1-P2) 


P2(8-8p2+3pS 

"  ud-02)5/2 


Hence  expanding  about 


p  we  get 


(3.1.33)  f(r)  = 


(1-P2): 


1  + 


p(2-p‘ 

2( 1-P‘ 


(r-p) 


p2(8-8p2+ 


8( 1-P‘ 


^2  ^  (r_p)2  +  o(n  ^//2) 


3/2  < 


Consider 


<p(r)  =  (1  +  p2  -  2pr )  a  . 


Differentiating  with  respect  to  r  twice  we  get 


> 


cp'(r)  =  ( -a)  ( l+p2-2pr )  a  ^(-2p)  =  2pa(  l+p2-2pr )  a  1 


and 
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cp"(r)  =  4p2a(a+l)  (l+p2-2pr) 


•a-2 


giving 


cp(p)  =  (l-p2)“a  ,  q)'(p)  =  2Pa(  1-P2) 


-a- 1 


and 


<PM(P)  =  Up2a(afl)  (1-P2) 


■a-2 


Thus 


(3. 1.3*0  <P(r)  : - - — T 

(1-P2)a 


1  +  (r‘p)  +  (r_p)2  +  0(n‘3/2> 


From  (3.1.33)  and  (3*1- 3*+)  we  get 


(3.1.35)  f(r)  <p(r)  = 


1 

H 

r  Pi2-P2) 

2pa 

P2)^ 

2(1-P2) 

1-P 

+  j 

p2(8-8p2+' 

1°  ' 

+i 

8(1-P2); 

+ 

O 

3 

l 

i\y 

J 

9 

The  coefficient  of  (r-p)  will  be  zero  if 


p(2-P2)  2pq 
2(l-pa)  +  1-P2 


=  0  , 


giving 


a  = 


P2  -  2 
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Substituting  this  value  of  a  in  (3.1*35)>  we  get 


f  (r)  <p(r)  =  - “ — ^  |"l  +  q(  j8-8p2+ 3P4+8( 2-P2) (^^)+l6(^j~) ( 

/t  L 


P2+2' 


(1-P2) 


+  0(n"3/2) 


(1-P2) 


p5  [x  +  8(  i-p^'l'g  {^p^p^-s+Sp^pN-p4-!* 


+  0(n"3/2) 


(1-P2) 


S  [l  +  ^'l-p'afg  (r'c)S  +  °(n*3/,2)l 


(3.1.36)  .*.  f(r)  cp(r)  = - ~ — ~ 

(1-P2)4 


p2(  p^-2) 

U(1-p2) 


(r-p)2 

1-P2 


+  0(n  3/2) 


Now  let 


(3.1.37)  x(0 


1  -  r2 

1  -  2pr  +  p2 


in  x(r)  =  in  (l-r2)  -  in(l-2pr+p2) 


Differentiating  with  respect  to  r,  we  get 


*'(r)  _  z£L 
x(r)  '  l-r 


2  + 


2p 

l-2pr+p2 


(3.1.38) 


■ 
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Differentiating  with  respect  to  r  once  again,  we  get 


(3.1.39) 


x(0  x"(0  -  [x'(r)]2  (l-r2)(-2)-(-2r)(-2r)  2p  2o  __  . 

- - (l-r2)2  +  (l-2pr+p2)2 


From  (3.I.37),  (3.1.38)  and  (3. 1.39)  we  get 

x(p)  =  1.  x'(p)  =  0  and  X"(p)  =  " 


Hence 


X(r)  =  1  '  Y~2  '  +  0(n'3/2)  . 


{ 


P2(2-P  ) 
1_r2  1  U(l-P2) 

l-2pr+p2| 


_  _  P5(2-pS 

- 1  mtT5) 


+  0(n'3/2) 


U-P2) 


Using  (3.1.36)  we  get 


f (r)  cp(r)  = 


(1-P2) 


2  1  l-2pr+p^ 


P2(  2-P  v) 
l_r2  7  ^(1-P2) 


{l  +  0(n-5A)l  . 


Hence  (3. 1.29)  gives  us 


a-i 


P2-2 


(3.1.U0)  h(r)  ^  K'  — D-r2)  ( l+P2-2pr )  h  j  i.3or+p? 


P2(2-P  ) 
1_r2  X  ^(l-P2) 


( l-2pr+p2) 


a- 1 

2 


} 


*  -|l  +  0(n"5//2)  j- 


K'  ( l-r2)‘ 


n  P2(  2-P 
~ " 1  +  1/1  ri? 


U( l-p‘ 


} 


n  _  .  / P2-2\  P2(  2-P  ) 

(l-2pr+p2)2  1  +  Wl-o'V 


1  +  0( 


I'5/2)}  , 
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where,  again,  K'  is  an  adjusted  normalizing  constant. 


Now  let 


(3.1.41) 


N  = 


P2(2-P4 
n  +  2(l-P2 


Hence 


—  -  1 


(3. 1.42) 


h(r) 


K1  ( 1-r2) 


fel  {  1  *  0("'3/a)  }  • 


(l-2pr+p2) 


2  4  2 


To  renormalize  h(r),  we  consider 


I  = 


I-  1 

(1 -r£)2 _ 


( l+P2-2pr) 


Nip 

2  "  2  "  4 


5-  dr 


Let  u  = 


1  +  r 

2 


Then 


and 


1  - 

(2u  - 

l)2 

1  - 

(4u2 

-  4u  + 

1) 

4u(  1 

-  u) 

9 

1  + 

P2  - 

2p(2u 

-  1) 

(l  + 

P)2 

-  4pu 

(1  + 

P)2 

{■- 

4p 

(1+p)2 

ul  • 


■ 
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Hence 


.  r1 


J 


N  .  N  .  „  .  p‘ 

2N_1  u2  (1-u)2  (1+p)  +  +2  -{  1 


{ 


M. 


(i+p); 


•i 


N  1  p* 
2  +  2  +  4 


du 


N  - 1  ..  N 

2  (1+p) 


M  i  P  ,  N  . 
-N+i+2  pi  ^2-l 

*  J  n 


(1 


N  - 

■u)2  '  { 


Nip' 


1  - 


^-u}: 


2  +  2  +  4 


(1+p) 


du 


Thus  (using  [41],  page  293*  example  number  l),  we  get 


'2  ,Nv  n/N< 


„N-1 


I  =  2  "(l+p) 


- (n- i-—  )  r (§)  r© 


2y  ^2'  r  AN  1  p2  N  4p  \ 

r(N)  V  2  2  4*2’  *  (i+p)2y  ’ 


where  F(a,b;c;z)  is  the  usual  hypergeometric  function. 

Using  the  duplication  formula  for  the  gamma  function,  ([41],  page 
240),  we  have 


J  r(H)  =  s"-1  r(|)  r(Sii)  , 


giving 


r(|)  r(|)  /n  r(|) 


TW 


2n-i 


Hence, 


-(N-  1-  §  ) 
I  =  sT*  (1+p)  ^ 


r(“)  / 
_ SL  F(^n 

r(N±l)  U 


4p 


1  .  e2  s. .  „  • 

2  k  •  2  ’  ’  (1+d): 
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Therefore  the  renormalized  density  function  is  given  by 


(3.1.43)  h(r)  'v  K(N,p) 


_ Uzjdf _ 

N  _  1  _  p 
(l+p2-2pr)2  2 


where 


(3.1. Vt)  K(N,o) 


_ (1+P 

'Tic  T(|) 


N-l 


N 


r(^) 

P2  N  . 

h  ’  2  * 


N 


(up)2 ; 


and 


»  — “S# 


When  p  =  0  ,  we  have 


K(N,o) 


=  K(n,0) 


p=o 


_1_ 

\/rt 


r^) 

f(|) 


F(Sti  ■  |  ;  n  ;  °) 


n/tc 


r(^) 

r(|) 


and 


r(£±A) 

(3. 1.1*5)  h(r)~-i- - f— 

r(l) 


(1  -  r2)2  |l  +  0(n'5/2)|  . 
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For  p  =  0  ,  (2.2.42)  gives  us 


K(N,P) 


i  r(ST+  D 

=°  =  A  r("  -  1+  x) 


_i_ 

A 


r(ffi 

r(|) 


and  from  (2.2.41)  we  get  the  renormalized  density  function  of  Daniels'  as 


h(r)  ~  — 


r(“*> 


&  r<f> 


n^.2 

2 


(l-r2)  2  |l  +  0(n'3/2)| 


so  that 


(3. 1.^6)  h(r )  ~  ~- 

vn 


r(^) 

r(|) 


(l-r2)2  -|l  +  0(n"5/2)|  . 


The  identity  of  (3.1,45)  and  (3«l.ii6)  shows  that  our  result  is 


the  same  as  that  of  Daniels'  for  p  =  0  and  thus  provides  a  sort  of  check. 
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CHAPTER  IV 

DISTRIBUTION  OF  AN  ESTIMATE  OF  THE 
SERIAL  CORRELATION  COEFFICIENT 
UNKNOWN  MEAN 


We  consider  the  same  Pre-stationary  Linear  Markov  Process  as  in 
Chapter  III  and  discuss  the  case  of  the  unknown  mean  in  the  present  chapte 
We  estimate  the  unknown  mean  x  by 


(U.l.l) 


x  = 


‘5X1+X2+X3+* 


,+x 


n-1 


+-^x 


n 


n 


and  p  the  serial  correlation  coefficient  is  estimated  by 


(U.1.2) 


where 


J 


(U.1.3)  c  =  (x1-x)(x2-x)+(x2-x)(x  -x)+. . .+(xn_1-x)(xn-x) 
and 


(k.l.U)  Co  =  ■|(x1-x)2+(x2-x)2+.  .  .  +  (xn_1-x)2+i(xn-x)2 
Now 

C  =  (x1-x)(x2-x)+(x2-x)(x  -x)+. . .+(xn_1-x)(xn-x) 
=  x1x2  -  x(xt+x2)  +  X2 
+  xgx^  -  x(x2+x^)  +  X2 


+  X  .X  -  x(x  _+x  )  +  x£ 
n-1  n  '  n-1  n' 


'I 


■ 
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=  c  -  2x(ix1+  x„  +  . , .  +  x  +  ix  )  +  (n-l)x2 
c  1  2  n-1  ^  n 

-  c  -  2x(n-l)x  +  (n-l)x2 
or 

(4.1.5)  C  =  c  -  (n-l)x2 

and 

Co  =  i(xrx)2  +  (x2-i)2  +  ...  +  (*n.ri)a  +  #(xn-S)2 

=  ^x2  -  X^X  +  -|x2 

+  x2  -  2x^x  +  X2 

+  x2  -  2x,x  +  x2 

+ . 

+  x2  -  2x  ,x  +  x2 
n-1  n-1 

12  ”  1  “  p 

+  -kx  -  X  X  +  2X 

^  n  n  * 

=  c0"2x(|x1  +  x2  +  x3  +  •••  +  xn_i  +  4xn)  +  (n-l)x2 
=  c  -  2x(n-l)x  +  (n-l)x2 
or 

(4. 1.6)  Cq  =  co  -  (n-l)x2 


The  joint  moment -genera ting  function  of  C  and  Cq  is 


M(T  ,T)  =  E(TqCo+TC) 


n 


,00  T  C  +TC-ix,c”1x 
o  o  - — 


dx 


n 


-00  -00 


e 


dxldx2*  *  * 


■ 
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where 


1+p2 

-p 

0 

0 

0 

-p 

1+p2 

“*  p  «  •  • 

0 

0 

0 

-p 

1+p2  .  •  • 

0 

0 

0 

0 

0  a  a  a 

1+p2 

-p 

0 

0 

0 

-p 

1 

— 

Now 


T  C  +  TC  -  \  x1  C-1  x 


o  o 


-1 


=  To(cQ-(n-l)x2)  +  T(c-(n-l)x2)  x'  C  x 

=  T  c  +  Tc  -  \  x'  C"1  x  -  (n-l)x2  (T  +T) 
oo  *  —  —  —  /  \D/ 


(T  +T) 

=  -  i  x1  B  x  -  ° 


n-1 


x  m  m  x 


where 


B  = 


Up2-to  -(p+T) 
- ( p+T )  1+P2-2T 


0 


0 

0 


-(P+T) 


0 

“(P+T) 

1+P2-2T 


0 

0 


0 

0 


0 

0 

0 


(P+T) 


0 

0 

0 


l+p2-2TQ  -(p+t) 


1-T 


and 


"  (x1,x2,...,xn)  , 

Hi  =  iz  *  1  1  >  •  •  •  >  ^  t~2  )  • 


Or 


T  C  +TC--4  x'  C"1  x  =  -4  x' 
oo  c-  —  —  —  — 


B  + 


2(Tq+T) 


,  m  m 
n-1 - 


x 


Hence 


M(T  ,T)  =  (2n )"n//2  f  ...  f  e 


o°  -|  x' 


B  + 


2(Tq+T) 


- : —  m  m 

n-1 - 


x 


-oo 


-00 


dx. dx  . . .dx 
12  n 


Hence 

(^.1.7) 


Now 


Hence 

(U. 1.8) 
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2(T  +T) 

M(T  ,T)  =  |B  +  - — 

o  '  n-1 


m  m 


(2 » ^ • • • » 1 »2 ) 


2(Tq+T) 

n-1 


1 

4 

1 

2 

1 

2 

1 

2 

•  •  • 

1 

2 

1 

2 

1 

4 

1 

2 

1 

1 

1 

•  •  • 

1 

1 

1 

2 

1 

2 

1 

1 

1 

•  •  • 

1 

1 

1 

2 

1 

2 

1 

1 

1 

•  •  • 

1 

1 

1 

2 

• 

1 

4 

1 

2 

1 

2 

1 

2 

•  •  • 

1 

2 

1 

2 

1 

4 

m  in' 

= 

S 

2 

S 

s 

s 

• 

•  • 

s 

s 

s 

2s 

2s 

2s 

• 

•  * 

2s 

2s 

s 

2s 

2s 

2s 

• 

•  • 

2s 

2s 

s 

2s 

2s 

2s 

• 

•  • 

2s 

2s 

s 

2s 

2s 

2s 

• 

•  • 

2s 

2s 

£ 

0 

s 

s 

s 

• 

•  « 

s 

s 

T  +T 
o 

n-1 


where 


s 


CD|OJ  CD  CD  CD  CD  CD  I  OJ 


-  84  - 


From  (3.1.3)  » 


(4.1.9)  B  = 


where 

and 


P+T 

q 

0 

0 


q 

p 

q 

0 


0  0 

q  0 
p  q 
q  p 


0  0 

0  0 

0  0 


p  =  1  +  p2  -  2T 


q  =  -(p+T)  . 


0  0 
0  0 
0  0 


0  0 
0  0 
0  0 
0  0 


p 

q 

0 


q 

p 

q 


0 

0 

0 

0 


0 

q 

l-T 


Hence 


(4.1.10)  |B  + 


2(T  +T) 
o 

n-1 


■  m  rrf|  = 


p+vf 

q+s 

s 

•  •  • 

s 

s 

q+s 

p+2s 

q+2s 

2s 

2s 

s 

q+2s 

p+2s 

•  •  • 

2s 

2s 

s 

2s 

q+2s 

•  •  • 

2s 

2s 

s 

2s 

2s 

•  •  • 

q+2s 

2s 

s 

2s 

2s 

♦  •  • 

p+2s 

q+2s 

s 

2s 

2s 

•  •  • 

q+2s 

p+2s 

s 

2 

s 

s 

•  •  • 

s 

q+s 

s 

s 

s 


s 

s 

q+s 

1'T,+f 


This  determinant  is  of  the  same  form  as  given  in  Appendix  II 


if  we  take  e  =  p  +  T  and  f  =  1  -  T 

o  o 


ro|w 
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Hence  from  (II. 14)  we  have. 


2(T  +T) 

(4.1.11)  |B  +  ° 


n-1 


rara'I  =  |\  - 


2s 


f  P+l+2q-n(p+2q)  \ 
l  (p+2q)2  1 


n 


s(p+2T  )2  s(p+2T  -2)2 

+  - —  X  .  +  - - -  Y  . 

2(p+2q)2  n‘X  2(p+2q)2  n_1 


+ 


Mr1 


s(p+2To)(2-p-2To) 

(p+2q)2 


n- 1 

q 


We  observe  that: 

1)  of  equation  (II.6)  is  the  same  as  |b|  . 

2)  X^  ^  of  equation  (II. 12)  is  equal  to  (-q)n  ^  A^_^  where 

^  is  defined  as  in  equation  (3.1»8). 

3)  Yn  ^  of  equation  (II, 13)  is  equal  to  (-q)  ^  where 

A'  ,  is  obtained  from  A  ,  by  substituting  e(=p+T  ) 
n-1  n-1  o' 

for  f (=1-T  )  . 

o' 

4)  p  =  1+P2-2Tq  gives  us  b  (as  defined  in  equation  (3.1,5)) 
equal  to  -  ^  . 

q 

5)  q  =  -(p+T)  gives  us  a  (as  defined  in  equation  (5* 1.5)) 

equal  to  -  —  . 

q 

Thus  q  =  -  —  and  p  =  -qb  =  —  . 

cl  3. 

We  now  proceed  to  evaluate  X  - ,  Y  .  and  various  other 

n-1  n-1 

expressions  so  as  to  be  able  to  evaluate  the  determinant  in  equation  (4.1.10). 
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Since 


X 


n- 


1 


A 

n- 


1 


from  equations  (2.2.20)  and  (3.I.I6)  we  get 


x  =  A-  [  ( l~z2n~2)  f  ( 1+z2)  ( 1-u-pr  )-rz(  l-P2))  _z2/ 
n  ^  a11  ^  L  l+p2-2pr-2u 


zn~\l-z2) 


or 


(4.1.12)  x  (UPa-2pr-2u)n^ 

n  ( l-z2)( l-2rz+z2)n 


x-  [(l-z2n  2){(l+za)(l-u-pr)-rz(l-p2)  }-z2(l-z^u  4)(  l+p2-2pr-2u)  ] 


2n-4, 


and  using  equation  (3,1.13)  we  get  from  (II.I3), 


,  __±T  m 

“-1  "  a""1  L  1- 


l-(b+aT  )z  ,  (b+aT  -z)z  , 
'  o'  n-1  v  o'  -n+1 
- z  +  -  Z 


1-  z‘ 


Now  using  equations  (2.2.22)  and  (3,1.11)  we  get 


b+aT  _  1+z2  ( 1+z2) ( u+pr )-rz( 1+p2) 

z( l+p2-2pr-2u) 


o  z 


( 1+z2) ( l+p2-2pr-2u+u+pr)-rz( 1+p2) 

z( l+p2-2pr-2u) 


or 


(4.1.13)  b+aT  g  .at^)lltO.!-Pr-u).  ..-.vzll±pS) 

z(  l+p2-2pr-2u) 


o 
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Therefore 


l-(b+aTo)z 


1 


( 1+z2) ( l+p2-pr-u)-rz( 1+p2) 

l+p2-2pr-2u 


} 


l+p2-2pr-2u-( l+z2) ( l+p2-pr-u)+rz( 1+p2) 

l+p^-2pr-2u 


or 


(4.1.14)  l-(b+aTQ)z 


-u-pr+rz( l+p2)-( l+p2-pr-u)z2 

l+p2-2pr-2u 


and 

(4.1.15)  b+aTo-z 


( l+z2)( l+p2-pr-u)-rz( l+p2)-z2( l+p2-2pr-2u) 

z( l+p2-2pr-2u) 


which  gives 
(4.1.16)  (b+aTQ-z)z 


l+p2-pr-u-rz(  l+p2)+(u+pr)z2 

l+p2-2pr-2u 


Hence 


Y 

n- 


1 


-fu+pr-rz(l+p2)+(l+p2-pr-u)z2]z2n  l+p2-pr-u-rz( l+p2)+(u+pr )z2] 

( l+p2-2pr-2u) ( l-z2)zn  1  a11  ‘ 


r  pn  •  2 

-fu+pr-rz( l+p2)+( l+p2-pr-u)z2)z^  +  f l+p2-pr-u-rz( l+p2)+(u+pr )z2} 


( l+p2-2pr-2u)( l-z2)z 


n-1 


zn  1  ( l+pa-2pr-2u)n  1 


( l-2rz+z2) 


n-1 


or 


2n-2 


(4.1.17)  Yn_^  =  [- {u+pr-rz(l+p2)+(l+p2-pr-u)z2}z 


+  {l+p2-pr-u-rz( l+p2)+(u+pr)z2} ] ( l+p2-2pr-2u) 


n-2 


(l-2rz+z2)n  (1-z2) 


x 
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Now 


(4.1.18) 


p  +  2T  =  1  +  p2  -  2T  +  2T 
o  o  o 


=  1 


+  P£ 


> 


(4.1.19) 


p  +  2T  -  2  =  -1  +  P2  , 
o  * 


P  +  2q 


2 

a 


b  -  2 


a 


z(  l+p2-2pr-2u)  /  JL 
l-2rz+z2  'Z  +  z 


2) 


or 


(4.1.20) 


P  +  2q 


( l-z)2( l+p2-2pr-2u) 

l-2rz+z2 


and,  from  (2.2.21)  and  (2.2.22), 

T  +T 
o 


1  f  ( l+z2)(u+pr)-rz( 1+p2)  z( l+p2-2pr-2u)-p( l-2rz+z2) 
n-1  l-2rz+z2  +  l-2rz+z2 


giving 


(4.1.21)  s  = 


u-p+pr+( l-r+pg-2u-par)z  +  (u+pr-p)z2 

(n-l)( l-2rz+z2) 
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Hence 


(4.1.22) 


l-(n-l)(n+2q) 

(p+2q)2 


1 


(n~l) 

(\-Tr:ll%)  (l+e2-2pr-2u)2 


f l-2rz+z2-(n-l)( l-z)2( l+P2-2pr-2u) ] T l~2rz+z2] 

(l-z)^  ( l+p2-2pr-2u)2 


(4.1.23) 


(-l)n'1(p+2To)(p+2To-2)qn"1 

(p+2q)2 


=  ( - 1 )n~ 1( 1+Pg) ( - 1+P5)1- l)n~ 1zn~ 1(l+p2-2pr-2u)n~ 1 

( l-2rz+z2)n~ 1  (Y-'Irz+z^)2  ( l+P2-2pr-2u)2 

(l-p^)zn  ^( l+p2-2pr-2u)n  ^ 

(l-z)h  (l-2rz+z2)n~5 


(4.1.24) 


r  p+2T  -] 

o 

2 

( 1+p2)  ( l-2rz+z2) 

2(p+2q) 

. 

2( l-z)2( l+p2-2pr-2u) 

(l+P2!2  (l-2rz+z2)2 

4(l-z)\ l+P2-2pr-2u)2 


and 


(4.1.25) 


2  1  _  f  ~( 1~P2)( l-2rz+z2)  1  2 

2(p+2q)  J  ”  [  2( l-z)2( l+p2-2pr-2u)  J 


(l-P2)2  ( l-2rz+z2)2 

4(  l-z)^(  l+p2-2pr-2u)2 
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Substituting  from  equations  (3. 1.18),  (4.1.12),  (4.1.17),  (4.1.18), 
(4.1.19),  (4.1.20),  (4.1.21),  (4.1.22),  (4.1.23),  (4.1.24)  and  (4.1.25)  in 
equation  (4.1.11)  we  get  after  simplification, 


B  + 


2(T  +T) 
o 

n-1 


m  m 


( l+p2-2pr-2u)n  ^ 

( 1-z2) ( l-2rz+z2)n 


:{l+p2-pr-u-r( l+p2)z+(u+pr )z2  }{l-pr-u-r( l-p2)z+(u+pr-p2)z2  } 


*  |l+02-2Dr-2u  +  2  /u-°tPr±(  hi 


[■ 


(1  'O' 


where  terms  which  are  relatively  0(n  ^)  have  been  neglected. 
Now 

1+p2-2Pr-2u  +  2  \  u~p.i-9,r+( A-r.1 2U±Q r. ) £+_( HiP+P O z 2  \ 

I  (1-z)2  i 


( l-p)2( l-2rz+z2) 

(1-z)2 


Hence 


(4.1.26) 


2(To+T) 

n-1 


( l+Pg-2pr-2u)T1  J _ 

( l-z)2( l-z2)( l-2rz+z2)n  1 


* ( ( l+p2-pr~u)-r( l+p2)z+(u+pr )z2}( ( l-pr-u)-r( l-p2)z+(u+pr-p2)z2} ( 1-p)2 


' 
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Therefore 

2(T  +T) 

(4.1.27)  M(u-rT,T)  =  M(To>T)  =  |B  +  n°x  m  m'  | 


n-1 

(l-z)(l-z2)2  ( l-2rz+z2)  2 

ml 

(1  +p2-2pr-2u)  2 

1 

{ ( l+p2-pr-u)-r( l+p2)z+(u+pr )z2)+2{ ( l-pr-u)-r( l-p2)z+(u+pr-p2)z2}+2( 1-p ) 


Combining  equations  (2.2.25)  and  (4.1.27)  we  get 

1  ml 

(4.1.28)  M(u-rT,T)  ||  ^  tirALX-**).  il\^Lz±^zl - 

( l-p)(  l+p2-2pr-2u)  2"" 

x.  _ 1 _ 

{ ( l+p2-pr-u)-r( l+p2)z+(u+pr )z2}2{( l-pr-u)-r( l-p2)z+(u+pr-p2)z2}2 


Differentiating  equation  (4.1.28)  partially  with  respect  to  u  and  setting 
u  =  0  we  get 

h  {,‘(u-rT*T)  I}] 

u=0 


1  ml 

(l“z)(l-z2)2  (l-2rz+z2) 

1  -  P 


2 


■' 


' 
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-n+3 

p  _  i  _  l. 

(n-5)(  l+p2-2pr )  {l+p2-pr-r(  l+p2)z+prz  2]  ?{l-pr-r(  l-p2)z+p(r-p)z2}  2 


-n+5 


.1 

2 


+  ^(  l-z2)(  l+p2-2pr)  { l+p2-pr-r(  l+p2)z+prz2)  {l-pr-r( l-p2)z+p(r-p)z2) 


l 

.2*1  “2 


=S±1 

4^( ] -z2)( l+p2-2pr)  2  {l+p2-pr-r( l+p2)z+prz2)  2{l-pr-r( l-p2)z+p(r~p)z2 ) 

(4.1.29) 

,2 


n-_ 

2 


n(l-z)(l-z  ).J  l-2rz+_z^)  —  j-  ^+p2_pr_r(  i+p2)z+prz2}"^{l-pr-r(  l-p2)z+p(r-p)z2  }“  2 


(1-p)  (l+p2-2pr) 


/i .  1  +  i . . . . 1 

^  n  2n  }  l+p2-pr-r( 


( 1+p2) z+pr z2  +  T-  p  r - r ( 1 - p  2 ) z+p ( r - p )  z 


T5 


Hence 


(4.1.30)  h(r)  ^ 


n 


( l-p)(  l+p2-2pr) 


ml 

$(z)  (l-2rz+z2)  2  dz 


where 


(4.1.31)  $(z) 


(l-z)(l-z2) 


1 

2 


[l+p2-pr-r( l+p2)z+prz2}2{ l-pr-r( l-p2)z+p(r-p)z2}2 


[*• 


1  ,  d“Z2)(  l+p2-2or) 
n  +  2n 


_ _  f  1  _  ,  _ 1 _ 

\  1+p 2-pr-r ( 1+p  2)~z+p  r  z  2  1-pr-r ( l-p2)z+p(r- 


Therefore 


roM 
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(4.1.32)  $(r) 


5 

_ (9-r-)lLir-!.).f _ f 

{ l+p2-pr-r  ( l+p2)r+pr ,r2)^ {l-pr-r ( l-p2)r+p(r-p )r2)2 


\  m  1  +  ( l-.r2,).Cl±P.^."2Pr )  f 
n  +  2n [] 


1  _  + _  1 

l+p2-pr-r( l+p2)r+pr ,r2  +l-pr-r(l-p2)r+p(r-p)r2 


'1 


{ l-r2+p2( l-r2)-pr( 1-r2) }2{ l-pr-r2( 1-pr ) }2 


T  { l+0(n  1)) 


3 

>2 


(  l"r)  (  l~r2)\ - T-  (  l+0(n-1)) 


( 1-r2) ( l+p2-pr )2( 1-pr ) 2 


Cl-r).(l-r2)2 

( l+p2-pr)2( 1-pr)2 


(l+Ofn'1)} 


As  in  Chapter  II  Section  2,  we  substitute  in  equation  (2.2.35)  f°r 
z  on  the  path  of  integration  and  expand  $(z)  as  a  power  series  in  z  -  r 
so  that  (4.1.30)  becomes 


(U.1.53)  fi(r) - 


n-4 

2 


00 


(2k) 


2«( 1-p)  ( l+p2-2pr ) 


m3 

2  k=o 


(r)(-l)k(l-r2)k  r  2k/ , 

(2k):  -  J  “  > 

-1 


2.-3 

2  2 


dw 
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Now 


Hence 


(4.1.31*) 


Therefore 


Using  an 


(4.1.35) 


where  IC 


xk  = 


f  («2)k  a-*2) 

'-1 


2  2 


du) 


■  I 


i  k-i  a-5 
k  2  /,  \2  2  , 
u  (1-u)  du 


0 


r(2|5) 

r{n±S!^2) 


I  = 
o 


r(i)  r(S|i)  /«  r(£jp) 


rf2!2) 


r(*f) 


T  .  1.3.5.  ...  (gk-3)(2k-l) 

k  (n-2)(n)(n+2) . . . (n+2k-4)  o 


irgument  similar  to  the  one  used  for  obtaining  (3. 1.29),  we  get 

n- 


h(r)  -v,  K  - 


2 


{ 1+0 (n  d)} 


(l+p2-2pr)  ( l+p2-pr)2( 1-pr )2 


is  an  adjusted  normalizing  constant. 


Proceeding  as  in  Chapter  III  for  obtaining  (3.IA0),  we  get 
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(4.1.36)  h(r)  ~  K 


(l-r)(l-r2) 

( l+p2-2pr  ) 


nZ 

2 


n- 

2 


Pg-2 


(  Up2-2pr ) 


(1-P2) 


J-i 

[1+P 


-r 


Paf  g-o;  I 

4( 1-p2) 


•2pr  I 


( 1+0 (n  jj 


K 


(1-P2) 


(.i-rJIl-rg). 


0S(2-Pk) 
2  +  4(1-PS) 


n 


( l+p2-2pr) 


P2-2  P2(  2-p1' ) 

4  +  4(l-p2) 


_  1 

{ 1+0 (n  2  )]  . 


Hence 


(If. 1.J7)  h(r)».K' 


N  j) 
2  ’  2 


( l+P2-2pr) 


N  P‘ 

2  '  1  "  4 


{ 1+0 (n  d)} 


where 


K’  = 


K 


(1-P2) 


and 


m  P2( 2-P  ) 

’  +  2(l-p2) 


Consider 


n 


r  ( l~r )( l-r2) 


N-3 

2 


J 


dr 


-1 


( l+P2-2pr ) 


N-2  P‘ 
2  "  4 


cu 
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l+r 

Substituting  u  =  — we  get 


N-5  N-5  N-3  N^ 

j  =  f1  2(l-u)  2  2  u  2  2  2  (1-u)  2  2du 

n  J  N-2  P2 

^  n2  P  k 

>n'2'5  r  up 
(1+o)  [1-a^pu 


.N-l 


(1+P) 


N-2-|£ 


/ 


N-1 

2 


u 


(1-u) 


N-l 

2  r  kp 


./Hli  _ef\ 

1  2  k  ' 


du 


Using  ( [ kl ] ,  page  293 >  example  number  1 ),  we  get 


.N-l 


(k. 1.38)  Jn  = 


(i+o) 


p‘ 

N-2-| 


r^)  r(S|l) 

r(N) 


c 


kp 

(1+P)2 


Using  the  duplication  formula  for  the  gamma  function  ([kl],  page  2k0)  we 
have 


r(^) 

r(N) 


n/tc 


2N-i  r(|} 


Hence 


sTt 


it 


n 


( 1+P) 


N-2-| 


r(^) 

rd) 


2  N-l 


N-2  P  _ 

2  “  k  *  2 


; N ;  fWP 
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Therefore,  the  renormalized  density  function  is 


(4.1.39) 


h(r)-vK(N,p)  ...U-.r)(.l-r.2l.g. 


.2 

[ l+0(n  2)] 


( l+P2-2pr) 


N-2  Pc 
2  “  b 


where 


(4.1.40) 


K(N,o)  = 


( 1+p  ) 


pc 

N-2-| 


r(|) 


r  TntF.  F  Nil  .  N  .  _ 

'/*  r(  2  5  FV  2  4  ’  2  ’  N  ’  7uSp„ 


and  F(a,b;c;z)  is  the  usual  hypergeometric  function. 
When  p  =  0,  we  have 


K(N,o) 


p=o 


=  K  (n,o) 


1 


r(|) 

rn-1, 


A  r(f) 


F(2ts  •  ;  n  ;  0 


and 


i  T(-) 
_1  v2; 

^  r(S|i) 


(4. 1.4l)  h(r)  -v 


r(-) 

1  k2j 


r(Sfi) 


-I 

(l-r)(l-r2)  2  [ 1+0 (n  2)] 
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For  p  = 


so  that 


(4.1.42) 


0,  (2.2.43)  gives  us 


h(r )  -v  — 


r  (“*) 


S*  r(2ji)  (n-1+1) 


n-1 


-1 


.  1 


(l-r)(l-r2)  2  [ 1+0 (n  2)] 


JL  2 
n 


!1  r(n} 
2  y2J 


Sl2 

2 


(l-r)(l-r2)  -  [  1+0 (n  *  )  ] 


r(Y) 


h(r) 


1  '2 


r(5) 


r(^) 


( l-r )( 1-r2)  2  [ 1+0 (n  2)]  . 


The  identity  of  (4.1.41)  and  (4.1.42)  shows  that  our  result  is 


the  same  as  that  of  Daniels'  for  p  =  o  and  thus  provides  a  sort  of  check. 
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APPENDIX  I 

SADDLE  POINT  APPROXIMATION 


In  the  present  appendix  we  discuss  section  2  of  Daniels'  paper  [5]. 
The  distribution  of  statistics  of  the  form 

c 

r  =  c  > 
o 

where  c  is  non-negative  is  to  be  considered.  Let  c  ,  c  have  the  joint 
o  o 

probability  density  f(co>c).  We  wish  to  find  the  distribution  of  r. 

The  Jacobian  of  the  transformation 


is 


c  =  rc  , 
o 


c  =  c 
o  o 


^(c0>c) 

d(co,r) 


Thus  the  joint  probability  density  for  cq  and  r  is 


c  f(c  ,rc  ) 
o  o  o' 


The  density  for  the  distribution  of  r  can  be  obtained  as  a  marginal  density 
by  integrating  out  cq  to  give 


(1.1) 


h(r)  = 


r 


J 


f (c  ,rc  ) 
o  o 


dc 


o 
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Let  M(T^,T)  be  the  joint  moment  generating  function  of 


and  c .  Then 


T  c  +Tc 

M(T0,T)  =  £(e  0  °  ) 


.  r 


j 


00  pOO  T  c  +Tc 

/  e°°  f(c,c)dc  dc 

./  o  o 


-00  -00 


By  the  Fourier  inversion  formula  we  have 


i  n  r  “(t  c  +tc ) 

f(c  ,c)  =  r-  ~TTg  /  /  M(T  ,T)  e  °  0  dT  dT 
'  o’  7  ( 2n l )  JJ°  ° 


The  paths  of  integration  are  along  the  imaginary  axes  in  the  Tq  and  T 
planes  from  -  i  oo  to  +  i  oo  or  along  any  allowable  deformation  of  these 
paths.  That  is,  along  any  path,  a-ioo  to  a+ioo,  such  that  no 
singularities  of  M(Tq,T)  lie  on  the  new  path  of  integration  or  between 
it  and  the  imaginary  axis.  Since  c  =  rc^  ,  we  have 


£( 


l 

(2jti)  2 


p  p  ”C  (T  +rT ) 

/  /  M(T  ,T)  e  °  °  dT  dT 

JJ°  ° 


Taking  the  linear  transformation 

u  =  T  +rT,  T  =  T 
o 

with  Jacobian 

S(To,T) 

S(u,T)  =  1  ’ 


we  observe  that  this  equation  becomes 
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f(c  ,rc  ) 
'  o  o' 


(2*02  // 


-c  u 

M(u-rT,T)  e  ° 


du  dT  . 


Then 


c  u 

f(c  ,rc  )e  dc 
o  o  o 


fj  M(u-rT,T) 


■c  u 
o 


du  dT 


c  u 
o  , 
e  dc 


The  Fourier- inversion  formula  gives 


r 

j 


00 


o 


1 

2jti 


c  u 

e  °  dcQ  =  M(u-rT,T) 


so  that 


f 

j 


00 


o 


c  u 

f ( c  ,rc  )  e  ° 
v  o  o7 


M(u-rT,T)  dT  . 

J 


Differentiating  under  the  integral  sign,  where  permissible,  with  respect 
to  u,  we  obtain 


c  u 

f(c  ,rc  )c  e  dc 
o  o'  o  o 


1 

2ni 


Putting  u  =  o  in  this  equation,  we  have,  by  equation  (i.l) 


h(r)  m*i  f  h  Mu-rT>T)] 


u=o 


dT 
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which  is  equation  (2.5)  of  Daniels'  paper  [5]. 


If  we  wish  to  transform  T  to  some  other  variable  z  by 
T  =  T( z , u) ,  then 

M(u-rT,T)dT  =  M[u-rT( z , u)  ,T( z , u)  ]  — - dz  . 


For  simplicity  we  write 


M[u-rT(z,u) ,T(z,u) ] 


^T(z.u) 

dz 


M(u-rT,T) 


dT 

51 


dz 


Then  we  have 


where  integration  is  along  the  transformed  contour  in  the  z-plane.  By 
differentiating  with  respect  to  u  under  the  sign  of  integration  and 
setting  u  =  o,  we  obtain 


(1.2) 


h(r>  m*i  f  b  [M<“-rT-T>  111 


dz 


u=o 


This  equation  known  as  Cramer-Geary  inversion  formula,  corresponds  to 
equation  (2.7)  of  Daniels'  paper  [5]  and  is  the  one  used  in  Chapters  II, 
III  and  IV. 


The  main  problem  of  this  paper  is  to  evaluate  approximately 
integrals  of  the  above  type  when  the  statistics  c  and  cq  are  calculated 
from  moderately  large  samples.  The  cases  considered  are  found  to  have 
integrands  which  can  be  written  in  the  form 


.  8b 

lo  8  - 


cp(z)  W(z)]n 


where  n  is  the  sample  size. 


z 


In  the  applications  considered,  cp(z)  may  be  expanded  about 
A 

z  where 


^'(z)  =  0 


and  the  resulting  series  integrated  term  by  term.  It  is  found  that  this 
leads  to  an  asymptotic  expansion  in  powers  of  (n  the  dominant  term 
of  which  is  used  as  the  approximation. 
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APPENDIX  II 


EVALUATION  OF  A  DETERMINANT 


We  are  interested  in  finding  the  value  of  a  determinant  of  the 


form 


(IX.l)^n 


e+2 

q+s 

s 

•  •  •  s 

s 

2 

q+S 

p+2s 

q+2s 

...  2s 

2s 

s 

s 

q+2s 

p+2s 

...  2s 

2s 

s 

s 

2s 

q+2s 

...  2s 

2s 

s 

s 

2s 

2s 

. . .  q+2s 

2s 

s 

s 

2s 

2s 

. . .  p+2s 

q+2s 

s 

s 

2s 

2s 

. . .  q+2s 

p+2s 

q+s 

s 

2 

s 

s 

•  •  •  s 

q+s 

r  s 

f+I 

This  is  the  same  as  the 


(n+l)x(n+l ) 


determinant 


(II. 2) 


1 

1 

2 

s 

1 

1 

•  •  • 

1 

1 

1 

2 

£ 

0 

S+2 

q+s 

s 

•  •  • 

s 

s 

2 

0 

q+s 

p+2s 

q+2s 

0  9  9 

2s 

2s 

s 

0 

s 

q+2s 

p+2s 

0  0  0 

2s 

2s 

s 

0 

s 

2s 

q+2s 

9  9  9 

2s 

2s 

s 

0 

s 

2s 

2s 

. . .  q+2s 

2s 

s 

0 

s 

2s 

2s 

. . .  P+2s 

q+2s 

s 

0 

s 

2s 

2s 

. . .  q+2s 

p+2s 

q+s 

0 

£ 

2 

s 

s 

•  •  •  s 

q+s 

.  2 
f+s 

Let  R  (x)  stand  for  the  operation  of  writing  the  elements  of 
the  i*"*1  row  after  adding  x-times  the  corresponding  elements  of  the  jth  row. 
Then  performing  the  operations 


• 

no  - 


and 


we  get 


(II. 3) 


-  s)  »  ^^(“2s)»  ....  ^(n-l)  _eS  ^  *  ^ni(  ~^s) 

R(n+l)l<-S> 


n 


1 

1 

2 

1 

i 

•  •  « 

1 

i 

1 

2 

-s 

e 

q 

0 

•  •  • 

0 

0 

0 

• 

ro 

CD 

q 

p 

q 

•  •  f 

0 

0 

0 

-2s 

0 

q 

p 

•  •  • 

0 

0 

0 

-2s 

0 

0 

q 

»  •  » 

0 

0 

0 

-2s 

0 

0 

0 

9  •  • 

q 

0 

0 

-2s 

0 

0 

0 

p 

q 

0 

-2s 

0 

0 

0 

•  f  * 

q 

p 

q 

-s 

0 

0 

0 

♦  ft 

0 

q 

f 

0 

i 

i 

l 

«  t  f 

l 

l 

* 

-s 

e 

q 

0 

f  f  t 

0 

0 

0 

-2s 

q 

p 

q 

4  i  • 

0 

0 

0 

-2s 

0 

q 

p 

t  4  * 

0 

0 

0 

-2s 

0 

0 

q 

t  9  ? 

0 

0 

0 

-2s 

0 

0 

0 

•  9  9 

q 

0 

0 

-2s 

0 

0 

0 

t  •  • 

p 

q 

0 

-2s 

0 

0 

0 

4  4  4 

q 

p 

q 

-s 

0 

0 

0 

«  9  • 

0 

q 

f 

n+1 


ci+1 


Ill 


1 

1 

2 

l 

l 

•  •  • 

l 

l 

1 

2 

0 

e 

q 

0 

•  •  • 

0 

0 

0 

0 

q 

p 

q 

•  •  • 

0 

0 

0 

0 

0 

q 

p 

•  •  • 

0 

0 

0 

0 

0 

0 

q 

•  •  • 

0 

0 

0 

0 

0 

0 

0 

«  •  « 

q 

0 

0 

0 

0 

0 

0 

•  •  • 

p 

q 

0 

0 

0 

0 

0 

•  •  • 

q 

p 

q 

0 

0 

0 

0 

•  •  ♦ 

0 

q 

f 

_  n+1 


or 


(II. 4) 


&  =  -2s  P  .  +  Q 

n  n+ 1  n 


where 


(II. 5) 


n+l 


and 


(II. 6) 


n 


0 

1 

2 

i 

l 

•  •  • 

1 

i 

1 

2 

1 

2 

e 

q 

0 

•  •  • 

0 

0 

0 

1 

q 

p 

q 

9  9  9 

0 

0 

0 

1 

0 

q 

p 

9  9  9 

0 

0 

0 

1 

0 

0 

q 

9  9  9 

0 

0 

0 

1 

0 

0 

0 

9  9  9 

q 

0 

0 

1 

0 

0 

0 

9  9  9 

p 

q 

0 

1 

0 

0 

0 

9  9  9 

q 

p 

q 

1 

2 

0 

0 

0 

9  9  9 

0 

q 

f 

e 

q 

0 

9  9  9 

0 

0 

0 

q 

p 

q 

9  9  9 

0 

0 

0 

0 

q 

p 

9  9  9 

0 

0 

0 

0 

0 

q 

9  9  9 

0 

0 

0 

0 

0 

0 

9  9  9 

q 

0 

0 

0 

0 

0 

9  9  9 

p 

q 

0 

0 

0 

0 

9  9  9 

q 

p 

q 

L° 

0 

0 

9  9  9 

0 

q 

f 

n+l 


a 


*  »  • 

•  t 
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Multiplying  the  sum  of  rows  2  to  (n+l)  of  P  .  by  (-l)/(p+2q)  and 

n+ 1  ' 

adding  to  the  first  row,  we  obtain 


n+1 


-n+1 

p+2q 

i  _  (_e  +q) 
p+2q 

0 

0 

•  it 

0 

0 

1  _  .la+f 

2  P+2q 

1 

2 

e 

q 

0 

•  •  • 

0 

0 

0 

1 

q 

p 

q 

•  •  • 

0 

0 

0 

1 

0 

q 

p 

•  •  • 

0 

0 

0 

1 

0 

0 

q 

•  •  • 

0 

0 

0 

1 

0 

0 

0 

•  •  • 

q 

0 

0 

1 

0 

0 

0 

•  •  • 

p 

q 

0 

1 

0 

0 

0 

•  •  • 

q 

p 

q 

1 

2 

0 

0 

0 

•  •  • 

0 

q 

f 

n+1 


n+1 


Multiplying  the 

sum  of  columns 

2  to 

n+1 

of 

Pn+1 

by 

p+2q)  and  adding  to  the  first 

column, 

we  obtain 

e+f+2q-n(p+2q) 

l  (e+q) 

r\ 

0 

r\ 

0 

i  (q+f) 

(p+2q)2 

:  P+2q 

u 

•  •  0 

u 

2  P+2q 

1  _  Ig±al 

2  P+2q 

e 

q 

0 

9  0  9 

0 

0 

0 

0 

q 

p 

q 

9  9  0 

0 

0 

0 

0 

0 

q 

p 

9  9  9 

0 

0 

0 

0 

0 

0 

q 

9  9  9 

0 

0 

0 

0 

0 

0 

0 

9  9  9 

q 

0 

0 

0 

0 

0 

0 

9  0  9 

p 

q 

0 

0 

0 

0 

0 

9  9  9 

q 

p 

q 

1  .  (,q+f) 

2  P+2q 

0 

0 

0 

0  9  9 

0 

q 

f 

or 


n+1 


. 


- 
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<“-7>  Pn+1  -  ~^t2d 


n 


2e-p 

+  2(p+2q) 


..  /  ..  xn-1  2f-p 

n  +  ^  2(p+2q) 


N 


n 


where 


(II. 8)  M  = 
'  n 


and 


(II. 9)  Nn  = 


— 

p-2e 

0 

0 

P-2f 

2(p+2q) 

0 

•  •• 

0 

2(p+2q) 

q 

p 

q 

0  0  9 

0 

0 

0 

0 

q 

p 

9  9  9 

0 

0 

0 

0 

0 

q 

9  9  9 

0 

0 

0 

0 

0 

0 

9  9  9 

q 

0 

0 

0 

0 

0 

0  9  9 

p 

q 

0 

0 

0 

0 

9  9  0 

q 

p 

q 

0 

0 

0 

9  9  0 

0 

q 

f 

p-2e 

0 

0 

0 

0 

P-2f 

2(p+2q) 

9  9  0 

2(p+2q) 

e 

q 

0 

0  9  9 

0 

0 

0 

q 

p 

q 

9  9  9 

0 

0 

0 

0 

q 

p 

9  9  9 

0 

0 

0 

0 

0 

q 

0  9  0 

0 

0 

0 

0 

0 

0 

9  9  9 

q 

0 

0 

0 

0 

0 

9  9  9 

p 

q 

0 

0 

0 

0 

0  9  9 

q 

p 

q  Jr 

Now  expanding  in  terms  of  the  elements  of  the  first  row  we  get 


(II. 10) 


M 

n 


P-2e 

2(p+2q) 


X 

n- 


1 


+ 


(-i  r1 


-P’2f 

2(p+2q) 


n-1 

q 


and 


(II. 11) 


N 

n 


P-2e 

2(p+2q) 


n-1 

q  + 


(-ir1 


-£=gL- 

2(p+2q) 


-  • 

0 

-  m 


where 


(11.12)  X  = 
n-1 


and 


(n.15)  t  . 


Hence 


(II. 1U) 


P  q  0 

q  p  q 

0  q  p 


0 

0 

0 


e 

q 

0 

• 

0 

0 

0 


0  0  0 
0  0  0 
0  0  0 


0 

0 

• . «  p 

q 

0 

0 

0 

...  q 

p 

q 

0 

0 

...  0 

q 

f 

q 

0 

...  0 

0 

0 

p 

q 

...  0 

0 

0 

q 

p 

...  0 

0 

0 

0 

0 

...  p 

q 

0 

0 

0 

•  •  •  ^ 

p 

q 

0 

0 

...  0 

q 

p 

■[* 

-  2s 

f  e+f+2q- 

n(p+2q)  1 

l  (p+2q)2 

j 

n-1 


n-1 


n 


s(2e-,p,)|  s(2f,-.Elg  Y 

+  2(p+2q)2  n-1  +  2(p+2q)2  n-1 

/  ,vn-l  s(2e-p_)(2£-£)  n-1 

+  (  l>  (P+2q)2  q 


Now 


(11.15) 

(11.16) 


=  f  , 


X2  =  pf  -  q2 


and  in  general 
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X  _=pX  -  q2  X 
n-1  n-2  n-3 


or,  equivalently, 


(ii. 17) 


x  „  -  p  X  ,  +  q2  X  =  0 
n+2  n+1  n 


Taking  E  to  be  the  forward  difference  operator  of  the  calculus 
of  finite  differences  defined  as 


(11.18) 


E  Xj  "  Xj+1  * 


we  may  write  equation  (H.17)  as 


(E2  -  pE  +  q2)Xn  =  0 


Solving  this  difference  equation,  we  find  that  the  roots  of  the 
auxiliary  equation 


z2  -  pz  +  q2  =  0 


are 


u 


£± 


Vp^- W 


and 


v  = 


p  -  sTp^-kq2 

2 


We  observe  that  u  +  v  =  p  and  uv  =  q2.  Hence  the  general  solution  of 
equation  (II.I7)  is 


(11.19) 


__  .  n  .  n 

X  =  k.  u  4*  v 
n  1  2 


Applying  conditions  (11.15)  and  (II.16)  to  equation(H,l7)  we  have 


Xx  =  k^u  +  k^v  =  f 


and 


k^u2  +  k^v2  =  pf  -  q2 


■ 


; 
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Solving  these  equations  we  get 


k 


1 


(p.-y).f  -  a2 

u(u-v) 


and 

k  =  ( u~ E ) f  +  q2 

2  v(u-v) 


Substituting  these  values  in  ( 11,19)  we  get 


t  =  rf(p-v)frial  u""1  +  t(u-p)f^g.  v"-1!  . 

n  |_  (u-v)  (u-v) 


Hence 


(II. 20) 


n- 1 


I-Cp-v)f-q2)  un-2  l(u-£)f+a^l  v 

(u-v)  (u-v) 


-1  ■ 


We  observe  that  by  interchange  of  rows  and  columns  we  can  make 
Y  ,  take  the  form  of  X  ,  where  instead  of  f  we  will  have  e  .  Hence 


n-1 


n-1 


(11.21) 


Y  _  f(p-v)e-qf}un"2  ( (,u~-p)e+q2]v 

n-1  (u-v)  (u-v) 


n-2 


Again 


(11.22) 


Q  =  eX  -  q2X  _  • 

xn  n-1  n  n-2 


Using  (II. 14),  (II. 20),  (11.21)  and  (11.22)  we  get 


4  =  [ 1  - 


(eX  .  -  q*X  _) 
v  n-1  M  n-2' 


. 
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s ( 2e-p) 2 
+  2(p+2q)2 


2(p+2q) 


2 

■2 


Y 

n- 


1 


n-l  s(2e-p)(2f-p) 

(p+2q)2 


n-l 

q 


(-l)11'1  sfge-p)(2f-p)  n-l 
(p+2q)2  q 


+ 


g 

e  "  2Tp~+2qTg  (Me+f+2q-n(p+2q)3 


(2e-p)2) 


X 

n- 


1 


+ 


(gf-p) 

2(p+2q) 


2 

2 


Y 

n- 


1 


-  q 


2s 


(p+2q) 


2  [e+f+2q-n(p+2q) } 


n-2 


or 


f-ir1 

(p+2q)2 


s(2e-p)(2f-p)qn  1 


+ 


e(p+2q)2  -  |  {4e[e+f+2q-n(p+2q) ] 


(2e-p)2] 


f (p-v)f-q2}un  2  +  f (u-p)f+q2)vn  2 

u  -  v 


2f - 


p+2q 


?) 


2 


f (p-v)e-q2)un  2  +  f (u-p)e+q2}vn  P 

u  -  v 


(p+2q) 


(p+2q)2-2s{e+f+2q-n(p+2q) ) 


((p-v)£-q2}un"^+f(u-p)£+q2V 


n-3 


u  -  v 


■ 
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Hence 


(11.22) 


where  u 


ns 

2(p+2q) 


((p-v)f-q2)^-  (u-l^q2)+{(u-p)f+q2}^-  (v-l+q2) 


s(2e-p)(2f-p) 

(p+2q)^ 


n-1 

q 


+  (p+2q)2  |^e(p+2q)2-2se(e+f+2q)  +  |(2e-p)2 


f  (p-v)f-qg]un~2  +  f(u-p)f+q2lvn  2 

u  -  v 


,  1  (2£-P)2  ri(p-v)e-q2}un"2 

2  Vp+2q'  [  u 

-  (■“2q")'^  £( P+2q ) 2-2s  ( e+f+2q) 


+ 


f (u-p)e+q2]vn  2 

v 


~((p-v)f-q23un~3+((u-p)£+q2}vn~3 

u  -  V 


and  v  are  the  roots  of  the  equation 
x2  -  px  +  q2  =  0  . 
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From  (3 


h(r) 


where 

<p(r)  = 


Hence 

h(r)  'v 


which  gives  us 

h(p)  ~ 


Now  by  Stirling' 


APPENDIX  III 


ACCURACY  OF  APPROXIMATION  OF  h(r) 


l.2b)  we  know  that 


n 


2  s/n 


Iklfl 


n- * 
2 


( l+p2-2pr) 


n-2 

2 


r(|  -1) 

*(r) 


liifl  <p"(r)  o(n'8) 

gtrry  J(7T +  0(n  ’ 


O-rf)2 


( l+p2-pr)2( 1-pr)2 


(l+p2-2pr)(l-pr+|  ) 
( 1+p^-pr )  (l-pr ) 


-  k) 


n 


2  s/jt 


llirf) 


n-2 

2 


n-2 


( l+p2-2pr)  2  ( l+p2-pr)2( 1-pr )2 


-1) 

r(^i) 


1  _  - <p."1.0-  of n”2^" 
1  2(n-l)  <p(r  )  +  °(n  'J 


1  + 


1r(  l+p2-2pr)(l-pr+|  ) 


n  l  ( 1+p^-pr )( 1-pr) 


+  0(n"2) 


,  r(“  -1) 
n  1  v2  ' 

2  A  (i-pa)*  r(aji) 


-1, 


1  +  OCn-1) 


s  Formula  [ Ul ]  we  know 


T(x)  =  e~X  xX"2  (2«)2 


■jl  +  0(x  1)j-  . 
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Thus 


or 


Hence 


r(|  -1) 

rC2-1) 


-5+1 


n  .  2 


/n  1x2  2  (at)*  (i  +  o[(|  -  i)"1]} 


e  -  (|-D 


n-1 


e  2  t2!1)2"1  (2*)*  (1  +  oKSji)-1]} 


=  £  <rf)'z  (£i )2_1  (1  +  0(n'1)) 


l  --1 

h  ,n-2N2  ,n-2 


-1, 


=  e2  (*=*)*  irf)  2  (1  +  0(n  A)) 


n-1 


"2 


1_ 

‘2 


=  e2  [1  -  2  [1  -  Pf]  ‘  [1  +  0(n” X)  ] 


1  1 


/v  ec  e 


*  [1  -  JTJ]**  n'*  !i  -  i]‘*  [1  +  0(n_1)  ]  , 


(III.1) 


r(|  -  i) 
r(^) 


«v  n 


i 

’2 


(III. 2) 


h(r)  at  r  =  p  is  0(n2) 


To  establish  a  bound  on  the  order  of  magnitude  of  £[(r-e)sl  , 


consider 


.  181 


(III. 5) 


because 


giving 


and 


whence 


§- 1 

(l-rg)2 


(l+ps-2pr)2  ( l+p2-pr)2( 1-pr)2 


=  /  -A-r/.  .  I' 

{  l+p2-2pr  J 


f-1 


( l+p2-pr)2( 1-pr)2 


( l+p2-pr )2( 1-pr)2 


1  - 


Jr^oiL 

l+p2-2rp 


since  1  -  r2  =  l+p-2pr-p2+2pr-r‘ 


l+p2-2pr-(r-p); 


(l-|ot)2(l+p2-|o|) 


Jr^-l 

l+p2-2pr  f 


l-1 


pr  <  | pr |  =  | pj  | r|  <  | p| 


1-pr  >  1  -  | P | 

1  +  p2  -  pr  >  1  +  p2  -  |  p  | 

_ i _  *  1, 

1-pr  -  1  -  | P | 


elcu 
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and 


1  +  p^  -  pr  -  1  +  o  -  p 


Let 


g(r)  = 


n 


2  \ fn 


( 1-r2) 


B.  i 

2 


r( 


r1 


( l+p2-2pr)  ( l+p2-pr )2( 1-pr )2 


r(2-iJi) 


Then  using  (ill.l)  we  get 


(III. 4)  g(r)  <  0(n2) 


(1-M)2  (1+P2-|P|)2 


{■- 


r  -  p 

M\p\ 


-  -  1 

v2  -i  2 


since 


(1  +  | P |  )2  =  1  +  P2  +  2 1 P | 


>  1  +  P2  -  2pr 


giving 


(r-p)2  (r-p)2 

( 1+ | p | ) 2  l+p^-2pr 


so  that 


!  .  JS^IL  <  !  . 

l+ps-2pr  \l+|p|/ 


We  know  that 


(III. 5)  0  <  g[(r-p)2]  =  f  (r-p)2  h(r)  dr 

sj 

-1 


-1 


ClCVJ 
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Now 


1-P 

ri+|p|  5-i 

J  (l+|p| )2  Us(l-u2)  (1+|p|)  du 

-(UP) 

ur?r 


-!z£ 

l+lp 


(1  +  lpl)3  f 

-(UP) 

l+lpl 


n 


-  1 


u2(  1-u2)  du 


-1 


=  (i  +  |p|)3 


r(|)  r(|) 


(|) 


n  i  n 

o-t  ”o 


(ar)2 


e 


«3> 


2 

2 


5 


2 

e 
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Jl 

=  0(n  2)  . 


This  result  combined  with  (III. 5)  suggests  that  the  distribution 

__i 

of  r  is  concentrated  about  r  =  p  within  a  range  which  is  0(n  2)  on 
either  side  of  p  . 

Therefore  from  (III, 3) 


n 

2n2 


-1 


(1-r2) 


n 


-1 


(l+pa-2pr)2  ( l+p2-pr )2( 1-pr)2 


<  1 _ 


(i+p2-|p|  )2(i-|p|  )■ 


T  e 


<i- »*»«{*  -  fea 


(tES 


2  /  \4 

Xr.fi  .  ^  1  (  JLz£ 


1+1  Pi 


(1+P2-|P|)2  (1-|p|) 


the  logarithmic  expansion  being  valid  because  |r  -  p|  <  1  +  |p|  giving 


r  -  p 

J  +  |P|. 


2 


<  1  . 


roio 
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Hence 


-  -  1 

(1-r2)2 


--1 

(l+p2-2pr)2  ( l+p2-pr)2( 1-pr)2 


(i+o2-|p|)2(l-|c|)2 


Thus 


(III. 6) 


.  (  r-p 

g(r)  <  0(n2)  e  2  \1+lpl 


a 


Suppose  | r  -  p |  is  0(n  ) 


Then 


'I  -  i>  fe)2  is 


and 


will  tend  to  zero  exponentially  as  n  -k»  provided  a  >  -  1 

Thus  from  (ill, 6)  and  (3,1,29)  we  see  that  g(r)  (and  hence  h(r))  is 

exponentially  small  for  values  of  r  outside  some  range,  on  either  side 

of  p,  which  is  0(n  2).  We  can  conclude  that  over  the  effective  range 

1. 

of  r,  r-p  may  be  considered  to  be  of  0(n  2). 


